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We further develop an approach to identify the braiding statistics associated to a given fractional 
quantum Hall state through adiabatic transport of quasiparticles. This approach is based on the 
notion of adiabatic continuity between quantum Hall states on the torus and simple product states — 
or "patterns" — in the thin torus limit, together with a suitable coherent state ansatz for localized 
quasiholes that respects the modular invariance of the torus. We give a refined and unified account 
of the application of this method to the Laughlin and Moore-Read states, which may serve as 
a pedagogical introduction to the nuts and bolts of this technique. Our main result is that the 
approach is also applicable — without further assumptions — to more complicated non-Abelian states. 
We demonstrate this in great detail for the level fc = 3 Read-Rezayi state at filling factor v = 3/2. 
These results may serve as an independent check of other techniques, where the statistics are inferred 
from conformal block monodromies. Our approach has the benefit of giving rise to intuitive pictures 
representing the transformation of topological sectors during braiding, and allows for a self-consistent 
derivation of non-Abelian statistics without heavy mathematical machinery. 



I. INTRODUCTION 

The discovery of the fractional quantum Hall (FQH) 
effect [1] has demonstrated that under the right condi- 
tions, an interacting electron system may enter a state 
with topological quantum order [2]. Laughlin's semi- 
nal treatment [3^ of this new state of matter predicted 
that FQH systems should display a unique and rich phe- 
nomenology beyond the quantized Hall conductance that 
had led to its discovery. This includes the presence of 
robust gapless chiral excitations at the edge, as well as 
fractionally charged bulk excitations, which furthermore 
obey fractional statistics. These characteristics allow one 
to distinguish a great wealth of different classes of FQH 
states including, possibly, ones where the statistics of the 
quasiparticle-type excitations are non-Abelian 21 IS]- The 
latter might facilitate a particularly robust route to fault 
tolerant quantum computing [6l[7]. Experimental control 
of the quasiparticle excitations and their possible utiliza- 
tion for schemes of quantum computing remain among 
the foremost challenges of the field. On the theoretical 
side, things are largely under control due to a remark- 
able correspondence between FQH trial wave functions 
and conformal blocks in certain rational conformal field 
theories . This formal correspondence has given rise to 
powerful field theoretic mappings that allow one, among 
other things, to infer the quasiparticle braiding statis- 
tics of the state in question [51. This inference, however, 
remains without microscopic justification in most cases. 
Such a justification requires showing that the procedure 
followed agrees with the result of adiabatic transport of 
quasiparticles (as defined by trial wave functions), which 
ultimately defines the statistics. This remains challeng- 
ing in many cases of interest. For Laughlin quasiparticles, 
this program has been carried out early on by Arovas, 
Schrieffer, and Wilczek [5]. The non-Abelian case has 
proven to be a profound technical challenge. For the 
Moore- Read, or "Pfaffian", state, a proof has recently 
been put forth [TO] following a series of insightful papers 



[TTHI3] further developing the plasma mapping [3] (for 
p + ip wave paired superfluids, a proof was given earlier 
in Ref. 12). 

Prior to that, a number of non-rigorous techniques had 
been developed to independently confirm the conformal 
field theory (CFT) result for the braiding statistics of 
the Pfaffian state. These techniques have the additional 
merit of recasting the non-Abelian statistics of the Pfaf- 
fian state into a different language that makes no use 
of conformal field theory or related modular tensor cate- 
gories. Such alternative languages might be particularly 
desirable in the Pfaffian case, which is believed to be rel- 
evant to the experimentally observed plateau at filling 
factor = 5/2 [14]. The first such approach is based 
on the interpretation jT5l - [T9] of the Pfaffian state as a 
p + ip wave Bardeen-Cooper-Schrieffer (ECS) state of 
composite fermions [20j . The second approach employs a 
strategy that has been successfully applied to interacting 
many body-systems since Landau's concept of a Fermi 
liquid, but only recently to electrons in the fractional 
quantum Hall regime. This strategy is to view the com- 
plicated interacting many-body state of interest as the 
adiabatic descendant of a simple, non-interacting state. 
As demonstrated in a series of recent works pTH55] . an 
adiabatic continuity with the desired features is given for 
a large class of FQH states by taking the thin torus or 
cylinder limit, which, as a formal limit, had been con- 
sidered earlier in Ref. [Ml This approach gives rise to a 
language of simple strings of integers, or patterns, that 
are associated to various incompressible FQH states and 
their quasiparticle-type excitations. The same patterns 
also play a central role in the recently discussed connec- 
tion between FQH states and Jack polynomials [55U37j . 
and are intimately related to the "patterns of zeros" de- 
scribing these states [55HlTj . The adiabatic continuity 
between such "thin torus patterns" and FQH states on 
general tori has been utilized in Refs. [23landl25lto derive 
the statistics of various Abelian FQH states and of the 
Moore-Read state, respectively. 
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For the Pfaffian state, there are thus a number of al- 
ternatives to the standard CFT method of obtaining the 
statistics, including rigorous results. However, the gen- 
eralization of these alternative methods to more compli- 
cated non-Abelian states has thus far been limited. It has 
been argued in Ref. [25] that the method using thin torus 
patterns and adiabatic continuity should in principle be 
generalizable to other non-Abelian states. At the same 
time, it is by no means obvious that the approach chosen 
there will always give rise to a sufficiently constraining 
set of equations to determine the statistics. That this is 
the case might naively be expected from the fact that the 
set of inequivalent solutions obtained in Ref. |^is identi- 
cal to those obtained |42l|43] from an assumed knowledge 
of the underlying CFT fusion rules, together with Moore- 
Seiberg polynomial equations [H]. Indeed, in both cases 
one obtains eight distinct solutions that are all related by 
overall Abelian phase factors. Furthermore, it is known 
that the thin torus patterns efficiently encode informa- 
tion about fusion rules [351 [3S] ■ One might thus conjec- 
ture that both methods generally produce the same re- 
sults. We will show below that this is not the case. To see 
why this need not be surprising, we note up front a num- 
ber of important differences between these two methods 
at both the conceptual and technical levels. A key step 
in relating fusion rules to braid matrices is to impose the 
validity of the "hexagon" and "pentagon" equations as 
they appear naturally in rational CFTs [H] . These equa- 
tions will not be explicitly enforced in our approach. In- 
deed, our framework requires no a priori assumption that 
the result of adiabatic transport of quasiparticles along 
braiding paths is purely topological in nature. Rather, 
this fact emerges naturally — together with the proper 
non-topological (Aharonov-Bohm) contributions — as a 
result of adiabatic transport. The assumptions under- 
lying these two methods are thus quite different. It is 
hence not immediately clear whether the "thin torus ap- 
proach" can be generalized to more complicated non- 
Abelian states. The main purpose of this paper is to give 
an affirmative answer to this question, and to flesh out 
this scheme of attack in greater generality, by applying 
it to the level fc = 3 Read-Rezayi state [47) . 

Even more generally, an efficient route from thin torus 
patterns to the statistics of the underlying state might 
be to use the information encoded in the patterns to 
relate them to some rational CFT. As a general dis- 
claimer, this is not what we will attempt to do in this 
work. Rather, we will use adiabatic continuity and re- 
lated assumptions to work out the statistics of a given 
quantum Hall state within a consistent framework that 
is independent of the field-theoretic assumptions that are 
traditionally employed in this field. The fact that we ob- 
tain results that are consistent with the CFT framework 
can thus be regarded as a non-rigorous, but independent, 
confirmation of the latter. We note that the connection 
between CFT and the patterns of zeros of a quantum 
Hall state, together with its implications about braiding 
statistics, have been studied in detail by Lu et al [48j. 



We believe that an added benefit to the method de- 
veloped here lies in the fact that an intuitive language 
is provided to describe non-Abelian statistics, and this 
language does not require the reader to have much back- 
ground in mathematical physics. To proceed, we thus 
give a refined and more detailed account of simpler cases 
already studied by this method. The basic ideas and un- 
derlying assumptions are fleshed out in Sec. [IT] where 
the simplest Abelian state, the v = \/2 Laughlin state, 
is studied. The generalization to the non-Abelian case is 
studied in Sec. [III[ where the Pfaffian state is considered. 
This section gives a rather more detailed and somewhat 
improved account of results first presented in Ref. US] 
In Sec. |IV[ we then show that exactly the same set of 
assumptions that sufficed to treat the Pfaffian case can 
also be used to determine the statistics of the j/ = 3/2 
(level 3) Read-Rezayi state essentially uniquely (up to 
an Abelian phase and complex conjugation). Our re- 
sults concerning the statistics of this state, as represented 
though thin torus patterns, are summarized in Sec. [V] 
The reader not interested in technical details, but rather 
more in this representation, is recommended to glimpse 
over Sees. |IB) |TT^ [HCj and |IIG[ to gather the nuts and 
bolts of the basic language used in this work, and then 
skip ahead to Sec. |v] (for the Read-Rezayi state) or Sees. 
HID HIE (for the Pfaffian state). In the remainder of 
the present section, we will review some basic formalism 
(I A), and then proceed to outline our general scheme of 



attack in |IB[ Some technical details are relegated to two 
appendices. 



A. Physics of LLL 

We work on a torus, identified as a rectangular 2D do- 
main of dimensions Lx and Ly subject to (magnetic) pe- 
riodic boundary conditions. We take the magnetic vector 
potential to be in Landau gauge, A = (0,x). The mag- 
netic length is set equal to 1 such that L^Ly — 2ttL, 
where L equals the number of magnetic flux quanta 
through the surface of the torus, which also equals the 
number of orbitals in the lowest Landau level (LLL). An 
infinite cylinder is obtained in the limit Lj; — >■ oo, with Ly 
kept finite. We first construct a basis of the LLL on such 
a cyHnder. It is given by iPn{z) = ^" exp(— — ^K^n^), 



where k = 2Tr /Ly, z = x -\- iy is the particle's complex 
coordinate, and ^ = exp(Kz). From the LLL states ip"^ on 
the infinite cylinder one can construct LLL states tpn that 
satisfy proper periodic magnetic boundary conditions (cf. 
Ref. I49p on a torus with finite L^ — kL. Fixing some 
unimportant overall phases, these boundary conditions 
read 



ipn{z + Lx) = e"^yipn{z) 
ipniz + iLy) = ipn{z) , 



(1) 



for the present gauge, and the orbitals fn{z) satisfying 
these conditions are then simply obtained by "repeating" 



FIG. 1. Landau level basis on the torus of dimensions x 
Ly. The orbitals ipn{z) form a ID periodic "lattice" in the x 
direction. Each orbital v?„(z) localizes a particle at a; = Kn 
while being delocalized in the y direction, leading to a "ring 
shape" geometry. Consecutive orbitals are separated by a 
distance 

the LLL orbitals of the cylinder along the x direction: 

Vn{z) ^Y^^Uod^) ■ (2) 
3 

For both the cylinder and the torus (with sufficiently 
large L^;), the n-th LLL orbital has the "ring shape" ge- 
ometry shown in Fig. [ij The orbital Lpn{z) localizes a 
particle in the x direction around x = nn to within one 
magnetic length, such that consecutive orbitals are sep- 
arated by a distance k. At the same time, each orbital is 
completely delocalized in y. We can view the orbitals 
as forming a ID periodic "lattice" along the x direction, 
with each orbital representing a lattice site. Note that 
we have Lpn+L{z) = ipn{z), and in this sense the "orbital 
lattice" satisfies ordinary periodic boundary conditions 
in n. A "thin torus limit" [1TH33J [SO] can be defined as 
K ^ 1. In this limit, the orbitals in the basis ([2| are well 
separated and have negligible overlap. 

It is clear that the choice of LLL orbital basis made 
above treats the x direction on the torus differently from 
the y direction. However, nothing prevents us from ex- 
changing the roles of x and y. A "dual" basis of states ^„ 
localized at y = nn (for k = 27^ /Lg;), encircling the torus 
in the x direction (Fig. [2| , can be obtained by formally 
"rotating" the basis, followed by a gauge transforma- 
tion, via 

Vn{z) = exp(ia;y)<^„(-iz)|«,_>j5 . (3) 

Alternatively, it can be shown (via Poisson resummation) 
that the Tp^ basis thus defined is related to the original 
basis ([2]) through a discrete Fourier transform, i.e 

Vn{z) = ^'YQyi^{-iY^^')Vn-[z) . (4) 
n' 

In the presence of the magnetic field, the single-particle 
Hamiltonian commutes with two magnetic translation 
operators, whose form in the chosen gauge is given by 

I — g-K(9x-iy) 

(5) 

f — p-I^Oy 



FIG. 2. The two Landau level bases (pn{z) and '^„{z) on the 
torus, shown here as a rectangular strip, a) The "original" 
basis, in which orbitals <p„{z) form a lattice in x with lattice 
spacing k — 2-k jLy and encircle the torus in y. b) The dual 
basis, in which orbitals ^„{z) form a lattice in y with lattice 
spacing k = 2-k jL^ and encircle the torus in x. 

The orbital bases and have simple transformation 
properties under the action of these two non-commuting 
translation operators. One easily verifies that 

ix</5n(^) = Vri^\{z) t^^niz) = e^"^„(z) (6a) 

tyipniz) = e"^"(^„(z) ty^„(z) = ^„+i(2) . (6b) 

All orbitals are thus invariant under the action of the 
operators t^^ and ty'", which represent magnetic trans- 
lations by Lx and Ly in the respective direction. This is 
equivalent to the observation that both the ipn as well as 
the ^„ orbitals satisfy the same periodic magnetic bound- 
ary conditions ([T]) appropriate to the gauge A = {0,x). 

We finally mention some other important symmetries 
of the problem under consideration. Inversion symmetry 
acts on wave functions via Iipiz) = ip{—z)^ and on the 
basis states defined above via 

Iipn{z) = f-n{z) , I'Pniz) = V-n{z) ■ (7) 

Similarly, while there is neither time reversal symmetry 
nor mirror symmetry in the presence of the constant mag- 
netic field, the combined symmetry does exist. We denote 
by T the antilinear operator that acts on wave functions 
via Tip{z) = 4>{—z*)*, and on basis states via 

T(p„{z) = (p-n{z) , T(^„(z) = 'Pniz) , (8) 

where the second equation follows from the first with 
Eq. Q. The reflectional part of r is obviously a reflection 
about the y axis. We can similarly define an antilinear 
operator f that performs a reflection about the x axis in 
conjunction with time reversal, and which acts on basis 
states via 

fipn{z) = tprriz) , f(^„(z) = ^_„(z) . (9) 

B. Outline of the method 

Our method of inferring the statistics of a quantum 
Hall state can be broken down into a few elementary 
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steps, which can in principle be apphed to any quantum 
Hall state that can be assigned well-defined ground-state 
patterns through a thin torus limit . Here 

we give a brief summary of the individual steps and the 
underlying ideas. A detailed development of these ideas 
will be given in the subsequent chapters, with applica- 
tions to Laughlin, Moore-Read, and fc = 3 Read-Rezayi 
states. 

1. Identify integer patterns characterizing the state 

It has recently become appreciated that a large class 
of trial wave functions can be characterized by simple 
sequences of integers. These patterns may be identified 
through various interrelated approaches. We we focus 
here on the approach based on the thin torus limit and 
adiabatic continuity pTH55] . This is not done for mere 
convenience, as it turns out that adiabatic continuity will 
play an essential role in the following. However, the pat- 
terns themselves can also be identified using various other 
methods, such as the Jack polynomial construction ^35r 
[32] and through "patterns of zeros" [38MT] . 

To be more specific, we will assume that the following 
program can be successfully carried out for the quantum 
Hall state in question. We assume that a Hamiltonian 
has been identified whose ground state lies within the 
desired phase. This Hamiltonian is assumed be local and 
to induce no Landau level mixing. The Hamiltonian can 
then be deformed into a limit that describes a thin torus, 
where either Ly <C 1 or Lj; <C 1. In this limit the ground 
state will approach a trivial product state of LLL orbitals, 
either in the ipn basis (for Ly <^ 1) or in the ^„ basis (for 
^ 1). These limiting states can be simply labeled 
by the pattern of occupancy numbers of successive LLL 
orbitals. Examples include 100100100 ... for the v = 1/3 
Laughhn state, or 202020 ... and 111111 .. . for the de- 
generate ground states of the Moore- Read state at v = 1. 
In these examples it has been demonstrated numerically 
[2n[22] that the deformation of the Hamiltonian into this 
limit can be done adiabatically, i.e., the gap of the incom- 
pressible fluid is maintained along the way. We believe 
that these observations can be extended, at the least, to 
all classes of trial wave functions for which local parent 
Hamiltonians, whose spectrum remains gapped in the 2D 
limit of an infinite plane, may be identified. 

2. Use adiabatic continuity to organize the space of 
elementary quasiparticle-type excitations 

It will further be assumed that adiabatic continuity, 
as described above, does not only hold in the sector 
spanned by the incompressible ground states, but also 
in the sectors obtained by adding quasiparticles or quasi- 
holes to the system. Specifically, in known cases of special 
trial wave functions for which parent Hamiltonians can 
be identified, the set of states obtained from the incom- 
pressible ground states by adding n > quasiholes may 
usually be characterized as the set of zero-energy states 
(or "zero-modes" ) of this Hamiltonian [51] . All examples 
discussed in the following will be of this kind. Adiabatic 
evolution from the 2D torus (by this we will always mean 
the regime L^, Ly ^ 1) into the thin torus limit will then 



take states with n quasiholes into states with n domain 
walls between different integer ground-state patterns [211 - 
I33| . For example, a state with n = 2 quasiholes in a i/ = 1 
Pfaffian may evolve into a thin torus state corresponding 
to the following occupancy pattern: 

2020202020202020111111111111111110202020202020202 . 

(10) 

A formal thin torus limit of the known zero-mode wave 
functions reveals the types of domain walls that represent 
quasihole states. The assumption of adiabatic continuity 
implies that at any aspect ratio of the torus, the zero- 
mode states will be in one-to-one correspondence with 
thin torus states of the form Eq. (10). One can thus 



define a complete basis of zero modes via adiabatic con- 
tinuation of the thin torus basis. The simple patterns 
that characterize the thin torus limit of a given basis 
state may still be used as "state labels" away from the 
thin torus limit. These labels carry important informa- 
tion about the transformation properties of basis states 
under magnetic translations. They also provide informa- 
tion about the change of topological sector for certain 
topologically nontrivial rearrangements of quasiparticles 
on the torus. The fact that these "thin torus labels" re- 
main meaningful, i.e., can be used to organize the space 
of zero-mode states away from this torus limit, will allow 
us to obtain information about the braiding statistics of 
the state, even though braiding statistics are well defined 
only on a (nonthin) 2D torus. 

3. Form coherent states describing localized quasiholes 
Individually, the adiabatically continued domain-wall 
states defining the quasihole basis described above do not 
correspond to states of well localized quasiholes. This is 
so because these states have a well defined momentum 
about the "quantization axis" , by which we mean the x 
axis when states are defined in the Ly ^ limit (where 
ipn orbitals are used to define product states), and the y 
axis when states are defined in the limit (where 

orbitals are used instead). Localized quasiholes ex- 
hibiting nontrivial braiding statistics are given by coher- 
ent state superpositions formed by states in the basis 
defined above. The general form of these coherent states 
is highly constrained by symmetries, non-commutative 
geometry (i.e., [x,y] cx i within a Landau level), and 
other consistency requirements that will be discussed as 
we go along. The validity of this general form may also 
be checked rather directly in the case of Laughlin states, 
see Sec. [TTJ 

4-. Determine transition functions describing a change 
of basis between dual coherent state descriptions 

A serious limitation of the coherent state ansatz men- 
tioned in the preceding step is that its validity is re- 
stricted to quasiholes that are well separated along the 
quantization axis (as defined above). However, there are 
two quantization axes at our disposal. These correspond 
to taking the opposite (mutually dual) thin torus limits, 
— >■ and Lx — ^ 0, respectively, giving rise to different 
ways of organizing the zero-mode space into adiabati- 
cally continued domain-wall states. These two ways of 
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organizing the zero modes are related by a modular S 
transformation of the torus, which essentially exchanges 
the roles of x and y. We will say that the corresponding 
different coherent state descriptions of local quasiholes 
are related by S duality. S duality allows us to write 
down a coherent state description for basically any lo- 
cal configuration of quasiholes. However, we will need to 
translate back and forth between mutually dual coherent 
state expressions along a braiding path. This change of 
basis is performed by matrix-valued transition functions, 
whose elements are sufficiently constrained by symmetry, 
topological considerations, and locality requirements to 
be discussed below. 

5. Adiabatically move the quasiholes along a braiding 
path 

The coherent state ansatz together with the transition 
functions allows the calculations of adiabatic transport 
of quasiholes along a given exchange path. In all cases 
studied, this confirms that the result of braiding is purely 
topological up to an Aharonov-Bohm phase, even though 
this is by no means a basic assumption made in our ap- 
proach. 

In the following, we will demonstrate the utility of this 
method for increasingly complex quantum Hall states. 
We start by discussing the simplest fractional quantum 
Hall state, the = 1/2 (bosonic) Laughlin state. 

II. THE LAUGHLIN STATE 

A. Thin torus limits 

Laughlin's v — 1/m wave functions [3] are the most 
elementary examples of a rich class of quantum Hall trial 
wave functions. These wave functions are generally char- 
acterized by a set of analytic requirements, the most ba- 
sic of which enforces that the wave function is entirely 
contained in the lowest Landau level (LLL). Laughlin's 
original construction of incompressible quantum liquids 
in a 2D planar geometry has been generalized by Hal- 
dane to states living on a sphere [52 enclosing monopole 
charges and to states on a torus 49J . The torus construc- 
tion has also revealed that the v — 1/m Laughlin state 
is m-fold degenerate on the torus, while it is nondegener- 
ate on the sphere. The nontrivial torus degeneracy was 
later understood to be the hallmark of topological order 
[2], and to be a necessary condition for the presence of 
anyonic excitations |53j . Here we focus on the torus. Let 
\'4>'^) , where c = . . . m — 1, denote the m incompressible 
Laughlin-type ground-state wave functions at filling fac- 
tor V = 1/m on the torus. We may expand the states 
in the basis of the LLL Fock space that is derived 
from the singleparticle basis (/3„: 

W) = ^ C{m^}\mi,m2---mL) . (11) 

{m„} 

Here, m„ denotes the number of particles in the state 
and we consider a system with a fixed number 



L = LxLy/2'K of flux quanta or LLL orbitals. For the 
time being, we will use Ly to parameterize the aspect ra- 
tio of the torus. The coefficients C{m„} depend on the y 
perimeter L y o f the torus. In the thin torus limit Ly — > 0, 
the states (ITT]) evolve into states dominated by a sin- 
gle pattern of occupancy numbers {tom}. E.g, the state 
with c — Q evolves into the Fock state |100 . . . 100 . . . ) 
(where dots indicate that I's are separated by m — 1 ze- 
ros), and states with c > are obtained by repeated ap- 
plication of the translation operator ■ is the many- 
particle version of the single particle translation operator 
tx discussed above, and acts on a thin torus pattern such 
as 100 . . . 100 ... as a right shift. For any value of the 
perimeter Ly, the Laughlin states are ground states 
of a "pseudopotential" Hamiltonian [551 [S3], whose ac- 
tion within the LLL explicitly depends on Ly. The evo- 
lution of the states with Ly can be understood as the 
adiabatic evolution of the ground states of the pseudopo- 
tential Hamiltonian H{Ly) as the parameter Ly is slowly 
changed. This has been studied in some detail for m = 3 
in Ref. [211 where is was shown numerically that the gap 
above the ground states never closes as a function of Ly. 

The thin torus states discussed here are formally iden- 
tical to the Tau-Thouless states proposed in Ref. [551 
When considered in the "2D-limit" L^ = Ly = oo, these 
states do not have long range charge density wave (CDW) 
order. In contrast, the thin torus states considered here 
can be characterized as ID CDW states breaking the 
translational symmetry of the system. This is so since 
in the thin torus limit, the LLL orbitals ipn are well sep- 
arated by a distance k = 2n/Ly (Fig. [2|, and the sym- 
metry breaking pattern of occupancy numbers becomes 
visible as a CDW modulation. The findings of Ref. JT] 
imply that the Laughlin states retain the CDW order of 
the thin torus limit on any torus with at least one of the 
dimensions L^, Ly finite. Related rigorous results have 
been discussed in Ref. [511 However, as long as both L^ 
and Ly are large compared to the magnetic length, the 
CDW order is exponentially small. The physics of the 
incompressible fiuid is thus quickly approached as L^, 
Ly become large, and in particular the notion of braid- 
ing statistics can be made arbitrarily well defined on a 
large but finite torus. This, together with the fact that 
the states on such a torus are adiabatically connected to 
simple product states sharing all their essential quantum 
numbers, is the foundation of the method discussed here. 

For simplicity we will now focus on the case m = 2, the 
bosonic 1/^1/2 Laughlin state with ground-state pat- 
terns 101010 . . . and 010101 . . . , respectively. The gen- 
eral case was worked out in Ref. 1231 However, here we 
will discuss an improved variant of the method, which 
was used in Ref. [25l to derive the statistics of the Pfaf- 
fian state. The two degenerate m — 2 Laughlin states on 
the torus are the unique zero-energy eigenstates of the Vq 
Haldane pseudopotential at filling factor v — 1/2. As in 
other cases where parent Hamiltonians for incompress- 
ible trial states are known, further zero-energy states ex- 
ist at smaller filling factors: The excitations associated 
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with elementary quasihole-type excitations are in one- 
to-one correspondence with the zero modes of the par- 
ent Hamiltonian at filhng factor v < 1/2. This is again 
true at any value of the perimeter Ly, and in particular 
the number of zero modes for any fixed number of con- 
stituent particles (electrons) N does not depend on Ly. 
We will extend the assumption of adiabatic continuity 
to the entire zero-mode sector. The thin torus limit of a 
Laughlin state with n quasiholes can easily be worked out 
directly from the Ly limit of the Hamiltonian [21], 
or from the same limit of the wave function on the torus 
or cylinder [SD] . A state with a single Laughlin quasihole 
evolves into a thin torus state that has a single domain 
wall between the two ground-state patterns. We can dis- 
tinguish domain- wall states in two "topological sectors" , 
according to the two possible phases of the charge density 
wave to the left and to the right of the domain wall, i.e., 

lOldOlOl... or OlOlOlOlO... . The ID domain walls can be 

I I 

ascribed a fractional charge by means of the usual "Su- 
Schrieffer" counting argument [S7]. This charge (here 
1/2) generally agrees [211 [29] with the charge of Laugh- 
lin quasiholes, as it should by adiabatic continuity. 

We introduce notation |a,c) for LLL product states 
with a domain wall at position a in topological sector c: 



\a,0) = 



\a,l) = 



loioioioidoioioioioio. 



.010101010100101010101 



(12a) 
(12b) 



The curved ket indicates that these are "bare" prod- 
uct states to be distinguished from states that have un- 
dergone adiabatic evolution, which we will discuss be- 
low. The number a is a half-odd integer labeling the 
domain-wall position relative to the LLL orbitals, such 
that a ± 1/2 are the orbital indices of the LLL orbitals 
adjacent to the domain wall. The two possible values of 
the topological sector label c distinguish the sequence of 
ground-state patterns in the two states of Eq. (13). It is 



worth noting that in principle, the topological sector is 
already determined by the value of exp(i7ra) and so the 
notation of Eq. (12 1 may seem slightly redundant. We 



find it advantageous, though, to include the topological 
sector information explicitly into the sector label, espe- 
cially with regard to more general cases discussed later. 

The above observations immediately generalize to 
states with two quasiholes, whose thin torus limits are 
given by product states corresponding to patterns with 
two domain walls. These states are labeled \ai,a2,c), 
with occupation number patterns for the values of c 
0, 1 given by: 



|ai,a2,0) = 
|ai, 02, 1) = 



.10100101010100101010. 

I I 

.OlOldOlOlOlOKiOlOlOl. 



(13a) 
(13b) 



We will always take ai to be less than 02, such that oi 
and 02 refer to the first and second domain wall, respec- 



positions are also subject to the constraint 

02 - ai = 1 mod 2 . (14) 

Again, the label c explicitly distinguishes the two pos- 
sible sequences of ground-state patterns, even though in 
principle this information is also contained in the values 
of exp(i7rai) or exp(z7ra2). The labels ai, 02, and c de- 
scribing a given two-domain-wall state are unique when 
the condition 



< oi < a2 < i 



(15) 



is imposed. Whenever the domain-wall positions sat- 
isfy (151, we will say that they are given "in the de- 



fault frame" . However, since we are working on the torus 
and LLL orbitals satisfy the periodic boundary condition 
ipn = ^n+Li it is desirable to admit domain- wall positions 
that refer to more general reference frames also. We thus 
define the states |ai,a2,c) for all ai, 02 satisfying 



ai < a2 < ai + L , 
together with the following identification: 
|ai,a2,c) = |a2 - L,ai,c') 



(16) 



(17) 



where c' = 1 + c mod m (here m=2). We will say that 
the domain- wall positions ai, 02 lie in an / frame if 



f<ai<a2<f + L 



(18) 



The standard frame is the frame. If necessary, repeated 
application of Eq. (17) allows one to transform domain- 



wall positions between different frames, where the roles 
of the first and second domain wall may be exchanged; 
whenever this happens, the topological sector label c 
changes also, as stated in Eq. (17). This fact follows 



from Eq. (14), since the value of c is determined by the 



value of the position of, say, the first domain wall modulo 
2, as discussed above. Note that L — 2N -I- 2 is even for 
states with two domain walls. The topological sector la- 
bel is therefore frame dependent. This is of a piece with 
the fact that the topological sector changes when one 
quasihole is transported around one of the "holes" of the 
torus, as we will discuss in detail below (see Fig.js]). The 
transformation properties of topological sectors under the 
exchange of two quasiholes along nontrivial loops (going 
once around the torus) arc thus encoded in the thin torus 
patterns. This is a key ingredient of the method pre- 
sented here, and sector transformation rules analogous 
to Eq. (17) will be of much importance especially in the 



lively. It is clear from Eq. ( 13 ) that the two domain- wall 



non-Abelian states to be discussed below. 



B. Delocalized quasihole states 

The notion of braiding is not well defined in the thin 
torus limit. In order for a well-defined statistics to emerge 
from an adiabatic exchange of quasiholes, throughout the 
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FIG. 3. Top: A possible arrangement of two domain walls 
is shown in a "repeated zone scheme" with the domain-wall 
positions marked by the bold 00 strings. The dotted line (red) 
marks the boundaries of the frame and the dot-dashed line 
(blue) marks a shifted frame, the 6 frame. Bottom: The 
second domain wall moves to a new position. When viewed 
from the frame, this domain wall moves across the frame 
boundary where it becomes the first domain wall in a different 
topological sector. Viewed from the 6 frame, the domain wall 
does not move across the boundary and the topological sector 
does not change. 



exchange the quasiholes must be spatially localized in 
both X and y, and at the same time must be kept away 
from each other at distances large compared to their in- 
dividual spatial extent. Both are not simultaneously pos- 
sible in the thin torus limit. Hence, in order to "braid" 
quasiholes through adiabatic transport, we will need to 
work with states that live not on a thin torus but on 
a full-sized torus with Lx,Ly both large. Formally, the 
assumption of adiabatic continuity means the following. 
There exists a family of unitary operators S{Ly, L'y) that 
describe the adiabatic evolution of the eigenstates (in 
particular the zero modes) of the pseudopotential Hamil- 
tonian at perimeter L'y into those at Ly. In particular, 

we define S{Ly) = S{Ly,0), the unitary operator that 
evolves thin torus states, Eqs. (12 1, (13), into states at 
finite L,.. We hence define the "dressed" or adiabati- 



cally evolved domain-wall states as the descendants of 
thin torus states via the operator S{Ly). In particular, 
for states with a single domain wall, we write 



\a,c,Ly) = S{Ly) |a, c) 



(19) 



where we will suppress the label Ly whenever no confu- 
sion can arise, using the regular ket to denote dressed 
states as opposed to bare domain-wall states. For suf- 
ficiently large Ly [and L^ = 2t: L/Ly), the states in 
Eq. ( 19 ) describe a quasihole immersed into a (here: 
h' = lj2) Laughlin liquid. The quasihole is localized in 
X around x = kg. However, it is entirely delocalized 
in the y direction. To see this, we consider the opera- 
tor Ty which is the many-body analogue of the single- 
particle translation operator ty discussed above. The 
bare domain-wall states are Ty eigenstates by construc- 
tion, with eigenvalues that are easily calculated from the 
pattern of occupation numbers. Since the pseudopoten- 
tial Hamiltonian commutes with the magnetic translation 
operators for any value of Ly, so does the adiabatic evolu- 
tion operator S{Ly). It follows that the dressed domain- 
wall states transform under magnetic translations in the 
same manner as the bare ones do. The states in Eq. ( 19 ) 



FIG. 4. A depiction of a two-quasihole dressed domain-wall 
state. These states are adiabatically evolved from the "bare" 
thin torus domain-wall states but live on the full-sized torus. 
The quasiholes described by this state are localized at some 
position in the x direction but delocalized in the y direction. 



those of their bare counterparts. It is clear that in such 
a state, the quasihole must be completely delocalized in 
the y direction (see Fig. |4|. Again, these observations 
can be extended to states with two quasiholes. 



\ai,a2,c,Ly) = S{Ly) |ai,a2,c) 



(20) 



via Aa; — kAo — 27r(a2 — ai)/Ly 
quasiholes in the state |ai, 02, c, Ly 



Here, two Laughlin quasiholes in the topological sector 
c are localized in x around xi ~ nai and X2 = Ka2, re- 
spectively, and are both delocalized in y. Note that the 
X separation between the two quasiholes depends on Ly 

The two delocalized 
will be uncorrelated 
as long as Ax is much larger than a magnetic length (set 
equal to 1). There are certainly no such correlations in 
the thin torus limit, and even at finite Ly both the corre- 
lation length of the incompressible fluid and the range of 
the interaction remain on the order of a magnetic length. 
As we increase Ly, the adiabatic evolution will therefore 
not induce any correlations between the two quasiholes 
as long as Ax ^ 1 remains satisfied. In this case, the lo- 
cal properties of each of the quasiholes will be the same 



as those of the single quasihole described by Eq. (191 



We emphasize once more that the adiabatically contin- 



ued domain-wall states in Eqs. (191 and (20 1 are neither 



simple product states, nor are they any longer "thin torus 
states" in any sense. Rather, the assumption of adiabatic 
continuity allows one to organize the zero-mode subspace 
into a basis labeled by ID patterns for any value of Ly. 
These patterns carry information about the properties 
under magnetic translations not only of the thin torus 
states, but also of their adiabatically descended counter- 
parts at finite Ly. Finally, it will be of some significance 
that, since the adiabatic evolution operator S{Ly ) is uni- 
tary, the dressed states of Eqs. (19) and (|20| are or- 



are thus still Ty eigenstates, with eigenvalues identical to 



thonormal, since the thin torus product states certainly 



C. Coherent states 

Individually, the members of the basis of zero-mode 
states defined above describe delocalized Laughlin quasi- 
holes. In order to analyze the braiding statistics of these 
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quasiholes, we need to form states where quasiholes are 
localized in both x and y. Laughlin has constructed an- 
alytic wave functions for such states [3], which are also 
zero-energy eigenstates of the pseudopotential Hamilto- 
nian. It must therefore be possible to write these lo- 
calized quasihole states as superpositions, or coherent 
states, in the zero-mode basis defined in the preceding 
section. 

We consider the single-quasihole case first. According 
to the above, it must be possible to write 



(21) 



for a state with a quasihole localized at complex coordi- 
nate h = hx + ihy. Here, we anticipate that to localize 
a quasihole, it is sufficient to include states of a single 
topological sector into the superposition, such that the 
localized quasihole state still carries a well-defined sector 
label. The left-hand side of Eq. ( [2T| ) is assumed to be 
a Laughlin single-hole state. Interestingly, as long as we 
assume that a zero-mode basis |a, c) with the properties 
claimed in the preceding section exists, the coefficients 
C{h,a) of this expansion are fully determined. To this 
end, we note that 



{a',c\a,c) = const x Sa,a' 



(22) 



The vanishing of Eq. ( 22 1 for a ^ a' follows since for 
different domain- wall positions the bare state |a',c) and 
the dressed state |a,c) have different Ty eigenvalues, as 
is easily seen by writing out the corresponding domain- 
wall patterns and calculating the action of Ty. On the 
other hand, the constant in Eq. (22 1 does not depend on 



a, since states with different domain-wall position a are 
related by repeated application of T^. From Eqs. (21) 
and (22 1, it follows that 



C{h, a) cx (a, c\tpc{h)) 



(23) 



We also expect only those states \a, c) to have any appre- 
ciable weight in the coherent state (21) whose domain- 



wall position X — na is close to the x position of the 
quasihole. We will assume that the coefficients C(/i, a) in 
this region are not affected by a change from periodic to 
open boundary conditions, as long as the torus is cut into 
a cylinder by a cut along y that is far away in x from the 
quasihole. In particular, it is clear from the discussion in 
Sec. II AI that such a cut would affect the local structure 
of the ipn LLL basis (in terms of which the states |a, c) 
have been defined) only by negligible amounts (for large 
Lx). For cylindrical topology, however, it is possible to 
evaluate the right-hand side of Eq. (23) explicitly. For 



definiteness, we explicitly write out the wave function 
for the Laughlin state \tpc{h)) on a cylinder of perimeter 



V'c(/i; zi...zn) = J|(^j - Tq) Wiit 



X e 



Here, = exp(KZi) and rj = exp^nh). Evaluating 



Eq. ( 23 ) amounts to evaluating the coefficients of "dom- 



inance patterns" in the polynomial of Eq. (24). This can 



be done using "squeezed lattice" methods discussed in 
Refs. [HIIIT]- This shows that the above wave function 
does indeed lie in a definite topological sector, as defined 
by the thin cylinder limit. [55] One finds: [55] 



\Mh)) = A/'^<?!'(/i,«;a)|a,c) 

a 

where 

(f){h, x) = exp 



^i{hy + n/K)x - ^{hx - xf 



(25) 



(26) 



and A/" is a normalization constant independent of h. The 
general form of the coherent state wave function Eq. ( 26 ) 



could have been guessed based on the following observa- 
tions. As a function of x, (j){h, x) can be interpreted as a 
"minimum uncertainty" coherent state of a particle con- 
fined to one spatial dimension. This is consistent with the 
fact that, after projection into a single Landau level, the 
X and y components of the position operator do not com- 
mute, but satisfy a position-momentum-type commuta- 
tion relation [x,y] cc i. y position can thus be regarded 
as X momentum, and vice versa. It is thus natural that 
the y position of the quasihole enters as a momentum- 
like phase twist in Eqs. (25), (26). On the other hand, 
as a function of h, 4>{h, x) looks like a lowest Landau 
level orbital of a charge 1/2 degree of freedom in the 
same magnetic field that is felt by the underlying elec- 
trons. These heuristic considerations will later allow us 



to generalize the coherent state form Eq. (25) to more 



complicated cases, where a direct derivation of the kind 
outlined here is not straightforward. 

The next logical step is to generalize the expression 
( [25] ) to states with two localized quasiholes. This is not 
difficult, as long as the two quasiholes at complex po- 
sitions hi and /12 are well separated along the x axis, 
i-e., h2,x — hi^x 3> 1. In this case, we can argue that 
the presence of the one quasihole does not influence the 
other, and the natural generalization of the coherent state 
Eq. ( 25 ) takes on the following form: 



\^c{hi,h2)) = A/"^ (j)(hi,Kai)cj){h2, ^02) |ai, 02, c) 



ai <a2 



(27) 

The function (j){h,x) is just as defined in Eq. (26). The 



prime in the above sum denotes the restriction of the 
domain-wall positions to values corresponding to the 
topological sector c. These are different for ai and 02, 



as a result of Eq. ( 14 ) . To be precise, we can define the 



topological sector c for two quasiholes via the following 
constraint on the domain- wall positions: 



ai = 2ni - 1/2 



02 = 2^2 + 1/2 + c (28) 



with integers 712 > ni. By default, the sum in Eq. (27) 



(24) is further restricted to domain-wall positions within the 
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The restriction to a different 
will be indicated by a subscript 



— hi^x ^ 1 holds, 
in a more formal 



(25) 



default frame, Eq. ( 
frame according to i 
/, \ipc{hi,h2)) f. 

For as long as the condition h2 x 
Eq. (27) can be inferred from Eq 
way, using assumptions about the action of local oper- 
ators on the adiabatically continued domain-wall basis. 
Locality arguments of this kind will play an important 
role in the following, and we will devote the next section 
to the development these arguments. 



D. Locality 

It is useful to formalize the assumptions that enter the 
factorized two-quasihole ansatz, Eq. (27). This naturally 



leads to general assumptions about the matrix elements 
of local operators within the zero-mode basis of adiabat- 
ically continued domain-wall states defined above, which 
will be of further relevance in much of the following. 
Let p{r) be a local operator, localized at some position 
r = {rx,ry). We will later consider p{r) to be the op- 
erator for the local charge density at r, but for now we 
wish to consider a generic (not necessarily single-particle) 
local operator. The action of this operator within the 
LLL Fock space depends on the aspect ratio of the torus. 
We first consider the action of p{f) on a bare domain- 



case, for the matrix element between these two states to 
be finite, it is clear that either Ui = hi for all i, or there is 
a single j such that aj ^ bj , with both naj and nhj in the 
vicinity of ■ Otherwise the patterns associated with the 
two states would differ even in orbitals that are far re- 
moved from Tx along the x axis, and their matrix element 
would be exponentially small. In particular, matrix ele- 
ments between states in different topological sectors are 
not possible (in the thermodynamic limit). Although at 
large Ly, the dressed domain- wall states |ai, a2, 03 . . . , c) 
are quite different from their bare counterparts, they still 
describe topological defects inserted into the torus at x 
positions nai. We will assume here and in the following 
that if the associated patterns of two dressed domain- 
wall states differ by many microscopic degrees of free- 
dom, then this is also true for dressed states themselves. 
In particular, if the patterns of two states differ in orbitals 
whose separation along the x axis is large compared to 
one magnetic length, we assume that their matrix ele- 
ment for any local operator will be negligible. For states 
with well separated domain walls, the observation made 
above for bare states then extends to their dressed coun- 
terparts. I.e., non-zero matrix elements are of the form 



p{ai,h) 



(29) 



where the ellipses represent other domain-wall positions, 
which must remain fixed but otherwise do not affect the 



wall state |ai, 02, as, . . . , c) (which for finite Ly is not an value of the matrix element, and again 



eigenstate of the pseudopotential Hamiltonian) . Quite 
obviously, the operator p{r) can only generate matrix el- 
ements between this state and some other domain-wall 
state 62, &3, . . . , c) if the associated pattern of orbital 
occupancy numbers differs only locally between these two 
states, for orbitals whose location lies within a magnetic 
length of rx ■ We will usually be interested in cases where 
the domain- wall positions Kai, ^02, was, . . . are all sep- 
arated by much more than a magnetic length. In this 



Tx to within a magnetic length. With these assump- 
tions, we can easily show that Eq. (27) describes two 
localized quasiholes, assuming that Eq. (25) describes 
a single localized quasihole. Let now p(r)he the lo- 
cal density operator. We consider the expectation value 
{il)c{hi,h2)\p{'r)\il^c{hi,h2)) for \h2x-rx\ > 1, and show 
that this expectation value reduces exactly to that of 
{il]c{hi)\p{r)\'4>c{hi)) , which we know to describe a single 
quasihole at position hi . Using Eq. ( 29 ) , we have 



{il}c{hi,h2)\p{r)\^c[hi,h2)) = A/""* ^ ^ (/)(ft.i, Kai)*0(/i2, Ka2)*(?!)(/ii, k5i)0(/i2, K62) (ai, 02, c| /5(r) |6i,62,c) 

ai,a2 61,62 

~ TV 2^ (j){hi,Kai)*4){h2, Ka2)* 4>{hi, Kbi)(j){h2, ^02) (ai, 02, c| p{f) 02, c) 

ai ,a2 ,61 

= TV (j){hi,Kai)*4){hi,Kbi)(j){h2,Ka2)*(j)ih2,Ka2) p{ai,bi) 

ai ,a2 ,61 

~TV^^ (p{hi, Kai)* (j){hi, Kbi) p{ai,bi) . 



(30) 



01,61 



In the above, the primes on the sums enforce all the nec- 
essary constraints such that the bras and kets correspond 
to domain-wall patterns in the topological sector c, cf. 
Eq. (28 ). In the second line, we have used that the matrix 



elements are diagonal in the second domain- wall position 



for \h2x — Txl ^ 1. Furthermore, for h2x — hix ^ 1 the 
constraint ai , 62 < 0,2 which the domain- wall positions 
obey becomes irrelevant due to the Gaussian nature of 
the (f) functions, and the sum over 02 in the third line 
simply yields the normalization of the single-quasihole 
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state, Eq. (25). The last line is, however, identical to 
{tpc{hi)\p{r)\iJc{hi)) ■ In words, this shows that when f 
is far away along the x axis from the second quasihole, the 
expectation value of p{r) reduces to that of a state with 
a single quasihole at hi. Similar arguments show that if 
f is far away along the x axis from the first quasihole, 
(p(r)) reduces to that of a state with a single quasihole 
at h2- Together, this shows that for h2x — ^ix ^ 1, the 



same in terms of the dual zero- mode basis Eq. (31 1 : 



state (27) describes two quasiholes localized at hi and 

h2. 



E. Dual description 



The coherent state expression (271 is in principle suited 



to calculate the Berry connection governing adiabatic 
transport [60H62] . However, as the arguments in the pre- 
ceding section have made clear, Eq. (27) can be expected 



to be accurate only in the limit of quasiholes that are 
well separated along the x axis. As can be seen in Fig. [7j 
the X separation of the quasiholes must vanish at some 
point for any exchange path, even though the absolute 
distances between the quasiholes remain large through- 
out. As a result, Eq. (27) is by itself not sufficient to 



fully calculate the result of adiabatic transport. 

The resolution to this problem lies in making use of 
the modular S invariance of the torus. Though we have 
so far only used the thin torus limit Ly — > 0, the physics 
must be invariant under an exchange of x and y. In doing 
so, we may now define a zero-mode basis by working from 
the limit — 0. In this limit, the zero modes of the 
pseudopotential Hamiltonian are domain- wall states that 
are occupation number eigenstates in the ^„ basis. The 
corresponding ground-state and domain-wall patterns are 
the same as those appearing in the Ly ^ limit, except 
that the associated charge density waves extend along 
the y direction of the torus. We denote the bar e domain- 
wall states in the ^„ basis with an overline, e.g. |ai, a2, c) 
for a two-domain-wall state. We now proceed in a man- 
ner that is completely analogous to the definition of the 
"original" zero-mode basis on a general torus, Eq. (20 1. 



To this end, we define a unitary operator S{Lx) that 
describes the adiabatic evolution of states from the "nar- 
row X limit" to a finite value of . We then define the 
general zero-mode basis for two-quasihole states via 



|ai,a2,c,Lj;) = S{Lx)\ai,a2,c) 



(31) 



where again, we will drop the label on the left-hand 
side whenever no confusion is possible. The states in 



Eq. (31) describe quasiholes that are localized in y but 



delocalized around the torus along x. Similar definitions 
are made for states with n quasiholes. We can form local- 
ized quasihole states in a manner completely analogous 
to Eq. (p7|). So long as Eq. (27) describes two localized 



quasiholes at positions hi and h2 for any aspect ratio of 
the torus, invariance of the physics under exchange of x 
and y implies that the following expression will do the 



4>{hi,Kai)4>{h2, Ka2)\ai,a2, 



ai <a2 



where 



(l>{h,y) = (j){-ih,y)\^. 



exp 



(32) 



;{hx+7:/K)y~ ^{hy - yf 
(33) 

hly > 1. 



2y ■ 



and Eq. ( 32 ) is now applicable to the case h 
We thus nave at least one valid coherent state expression 
for any configuration of the two quasiholes along the ex- 
change path shown in Fig. [7j At some points along the 
path, however, we will be forced to translate back and 
forth between the two coherent state expressions (27) and 
(32). This task is nontrivial. To see this, it is important 
to note that the topological sector label c has different 



meanings in the original zero- mode basis Eq. ( 20 ) and the 
dual zero-mode basis Eq. (31): in the former, it means 



that the state evolves into a well defined charge density 
wave product state in the limit Ly — ^ 0, characterized 
by a certain sequence of ground-state patterns separated 
by domain walls; in the latter, it means the same in the 
opposite thin torus limit, L^ — > 0. It will turn out that a 
state that carries a definite sector label c in the original 



basis, Eq. (20), is a superposition of states carrying dif- 
ferent topological sector labels in the dual basis Eq. (31 1, 
and vice versa. The same is true for the coherent state 
expressions Eqs. (27) and (|32|. While the relation be- 
tween the sets of states |'i/'c(^i, ^2)) and \ipc{h'i,h2)) is 
thus not diagonal in the topological sector label c, for 
given quasihole coordinates hi, /12 both sets span the 
same subspace, namely the space associated with having 
quasiholes localized at hi, h2- The relation between the 
states |V'c(^i, ^2)) and \4'c{hi,h2)) is thus diagonal in the 
quasihole positions, and we may write 



\ipcihi,h2)) 



E 



<,,(/ll,/l2)|V'c'(/il,/i2)) 



(34) 



Note that above, we had defined \ipc{hi, h2)) only for 
h2x > hix, and \tpc'ihi,h2)) only for h2y > hiy. While 
we will stick to these restrictions most of the time, 
we will generally let |V'c(^2, /ii)) = \'^c{hi,h2)) and 
\4'c'{h2,h\)) = |^/'c'(/ii, /12)) for convenience. This allows 
us to write relations such as Eq. ( 34 ) without distinguish- 



ing different cases. The transition functions u^^,{hi, h2) 
are then meaningful in regions where both \hix — h2x \ S> 1 
and \hiy — h2y\ 3> 1, since it is only in these regions 
where we have defined both \ipc{hi, /12)) and \'ipc'{hi, ft.2)) 
through coherent state expressions. The final technical 
obstacle is to sufficiently determine these transition func- 
tions from symmetries and topological considerations. 

To this end, we begin by distinguishing two regions 
of the 2-hole configuration space. Let a — sgnihix — 



h2x)sgn{h 



l2y) 



(T = ±1 then refers to first and sec- 



ond quasihole configuration in Fig. [5] respectively. We 
will first be interested in the "local" dependence of the 
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FIG. 5. The two possible configurations of two quasilioles, 
which are distinguished by the value of a = sgn(/!,i2, — 
h2x)sgn{hiy — h2y). Left: a — +. Right: a — 



transition functions on coordinates within each of these 
regions. Later we will use the fact that these regions 
are actually connected by "global" trajectories where one 
quasihole is taken around one of the holes of the torus 
(Fig.|6]). For now we will not allow these global moves. 
Within each of these regions, we now show that the lo- 
cal dependence of the u functions on coordinates is as 
follows, 



^{hij,hiy + h-2xh-2y) — 



where the parameters are complex constants and 
u{hi,h2) is the phase function ei'^'^^'^'^^y~^'^^^'^^y\ 

The hi,h2 dependence of m^^, can be locally deter- 
mined from the Berry connections. Using the coherent 
state expressions in Eqs. (27 1 and (32) on the full-sized 



torus {k, k <^ 1), the Berry connections can be calculated 
to be 



i ii^cihi, h2)\ V/ii_2 \^pc{hi,h2)) = -^(0, hi^^2x) 



(36) 



z('(/'c(/ll,/l2)|V/ii_J-0c(/il,/l2)) = ^{hly,2y,0)- 



An essential ingredient in the above is the fact that the 
zero-mode basis states we have defined are orthonormal, 
as explained at the end of Sec. |IIB| This is where the 
assumption of adiabatic continuity is crucial in our ap- 
proach. Obtaining Eq. (36) is then straightforward, since 



(35) 



in the limit k, k <C 1, the remaining sums can be replaced 
by Gaussian integrals. 

Let us consider an adiabatic process where one quasi- 
hole is fixed at hi and the other is dragged from /i2 to 
h'2 (which are both in the same region a). This process 
is described by a unitary operator, which acts separately 
on each term on both sides of Eq. ( 34 ) , yielding 



exp 




IV'c(/ii,/i2)> =X!"cc'('ii>'i2) exp [ z / dh- 




{hy,0) 



liPAhuK))- (37) 



The above equation may be compared to Eq. (|34| eval- 
uated at {hi, h'2) instead of (/ii,/i2)- This yields a rela- 
tionship between the u functions at these two locations. 



Ucc'i^i^K) = <c'('ii'^2) exp i / dh 



^^hhxhy 



= <^,(/ii,/i2) exp \^-i{h2^h'2y - h2xh2y) 

(38) 

where we used the fact tha,t{hy,0) — — {0, hx) + 'V hhxhy. 
In order to satisfy Eq. (38), the dependence of u on ft,2 
must be proportional to e^'^'^^^'^^y . Using a similar ar- 
gument in which the quasihole at ft,2 remains fixed while 
the quasihole at hi is moved, we find that the depen- 
dence of u on hi is proportional to e^'^'^'^'^'^'^y . Therefore 
the general form of the u functions is given by Eq. ( 35 ) . 



F. Symmetries and further simplifications 

With the above considerations, the transition functions 
w^^/ have been reduced to parameters ^JT^,, , of which there 
are eight at i/ — 1/2. We will now establish further re- 
lations between these parameters using symmetries and 



adiabatic transport along the "global" trajectories men- 
tioned above. 

First, we derive relations arising from properties under 
magnetic translations. The magnetic many-body transla- 
tion operators , Ty introduced above have the following 
effect on the dressed domain-wall states: 



\ai,a2,c) = \ai + l,a2 + 1, 1 - c) 
Tx\ai,a2,c) = e^^3"-'|ai,a2,c) 



Ty loi, 02, c) = e 



ai, 02, c) 



Ty\ai,a2,c) = \ai + 1, 03 + 1, 1 - c) 



(39) 



(40) 



where c — 0, 1, and rij is the orbital index of the orbital 
occupied by the j-th particle in the thin torus pattern as- 
sociated with the state. For the bare product states asso- 
ciated with these patterns, the above identities are direct 
consequences of Eqs. ^ for the single particle transla- 
tion operators. However, the properties under magnetic 
translations remain the same for the dressed states, as ex- 
plained in Sec. II B Note that the basis states |ai,a2,c) 
are eigenstates of Ty whereas T^, changes the topological 
sector label, and vice versa for the basis states |ai, 02, c). 
Equations ( 39 ) and ( 40 ) allow us to work out the prop- 



erties of the coherent states under magnetic translations. 
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The fact that both sides of Eq. (34 1 must transform 



the same way under these translations poses severe con- 
straints on the coefficients . Observing that for given 
domain- wall positions, 

I]?ij = ^i(^^ + c)-^(ai + a2), (41) 



it is a simple thing to verify the following properties of 
the coherent states under magnetic translations: 

\Mhi,h2)} = e-^«(''i«+''^«)+" |V;i_e(/ii + k>2 + k)) 

(42) 

Ty \Mhuh2)) = e*-^+'-- \^,{hl + IK, h2 + IK)) 

Ty\Mhi,h2)) = e3-(''-+''-)+-|^i_,(/ii + IK, h2 + ik)) 

(43) 

where we define A = (which in the present case, eval- 
uates to the integer L/2 = N + 1). 

We can use these translational properties to constrain 



the eight Ccc'^- recast Eq. (p54|) in matrix form 

|V'o(/il,/i2)> 



\tpo{hi,h2)) 
\Mhi,h2)) 



= u{hi,h2)'BF 



\MhiM)) 



, (44) 



where we have used Eq. ( 35 1 , and ii'' is the matrix with 
elements . Let us apply Ty to Eq. ( 44 1 . 



\MK.h'2)) 



u{h[,h'^)E'^{en<^. 



(45) 



The positions hj — hj + ik for j = 1, 2, and the u{hi, h2) 
function has been shifted by absorbing the spatially de- 
pendent phase in Eq. (43). If we compare Eq. (|45|) to 



Eq. (44) evaluated at the shifted positions (/i']^,/i^)7 we 
find that the two equations are consistent, provided that 
the 'EF matrix satisfies the following constraint: 



e (Tz^ 



(46) 



We can derive another constraint using the same logic 
after translating Eq. ( 44 ) with : 



(47) 



These two sets of equations constrain the "EF matrix to 
be of the following form. 



V2 I 



1 



-1 



(48) 



where is a pure phase, and the overall normalization 
factor 1/^/2 has been determined from the requirement 




C)h2b 

I 



FIG. 6. Different configurations o can be connected by drag- 
ging one quasihole along a global path. Initially, the two 
quasiholes at h\ and /12 are in configuration a = +. Keeping 
the quasihole at hi fixed, the quasihole at /i2 can be moved 
along one of two paths: path a, in which the quasihole at h2 
moves around the torus in the x direction to /i2a, or path 6, 
in which the quasihole at /12 moves around the torus in the y 
direction to /126. Both paths can be used to change the con- 
figuration a while keeping quasiholes well separated in both x 
and y. At the same time, the topological sector also changes. 



that S°' is a unitary matrix. Thus, after using transla- 
tions we have only two unknowns remaining, the overall 
phases and It is only the relative phase between 
the two that will have physical significance. 

In order to fix this relative phase, we will now drag one 
of the quasiholes in a two quasihole state along a "global 
path", i.e., a path where the quasihole disappears on one 
end of the standard frame (see Sec. |ll and Fig. [S]) and 
reappears at the other. The merit of such a path is that it 
connects the a — + and a = — configuration while main- 
taining both conditions \hix — h2x \ ^ 1, \hiy — h2y \ 3> 1. 



Let us consider the coherent state \'4>c{hi,h2)) , Eq. (27l, 
with two quasiholes in the topological sector c in the 
a = + configuration. We will drag the second quasihole 
along path "a" as shown in Fig. [6) We will do so by con- 
tinuously changing the position of this quasihole from a 
value hi, with well within the boundaries and L^, 



to a value with < < hix + L^. The default 
frame introduced in Sec. Ill Al is not suited to describe 
this process continuously. We thus choose an / frame as 
described in Sees. Ill Al and jll CI and consider the state 



|'(/'c(^i, ^2))/, i-e., the coherent state (27) with the sum 



restricted to the / frame. For this we choose a param- 



eter / such that Kf < hix < h\x < hi^^ < K{f + L) = 



kJ + Lx- Note that as long as the x position h2x of the 
second quasihole is well between hix and Lx, one has 
\4'c{hi,h2)) f = \ipcihi,h2)), where = denotes equality 
up to exponentially small terms. In this case the weight 
of both Gaussians in the coherent state is well contained 
within both frames, and so {ipcihi, ft.2))/ and |-0c(/ii, ft.2)) 
may be used interchangeably. However, as soon as h2x 
approaches Lx, we must work with |-0c(^i, /12)) /• In this 
regime, we will see that the coherent state \ipc{hi,h2)) f 
is identical up to a phase to the (default frame) state 
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|'0c'(^2 — Lx,hi)). That is, the second quasihole reap- 
pears on the left end of the standard frame, thus becom- 
ing the new 'first' quasihole (Figs. [3]and|6|. However, 
in the default frame the final state will be in a different 
topological sector with c' = 1 — c. At the same time, 
the quasiholes are now in the a = — configuration. This 
allows us to obtain one more relation between the tran- 
sition functions w^^, and their defining parameters . 

We first establish the precise relationship between 
\ipcihi,h2)) f and |'!/'c'(^2 — Lx^hi)), where exceeds 
by more than a magnetic length. One finds: 

\Mhi,h2)} f = e5^''^«^^+-''+-|V;i-e(/i2 - L,, h,)} j^l 

(49) 

where in the first identity we have passed to the f — L 
frame by straightforwardly plugging the identification 
(17) into the coherent state ((27|). The second identity 



follows from the fact that for h2x well exceeding L^, the 
states \i>c'{h2 - Lx,hi))f-L and |'(/'c'(/i2 - Lx.hi)) are 
again identical up to exponentially small terms, as dis- 
cussed above. 

Next we look at the comparatively trivial issue of how 
the dual state \^c'{hiTh2)) transforms along the same 
path, where h2 is again taken from ft,2 to Since the 
motion is chiefiy along the x direction, there is no need 
for a change of the frame for the |ai,a2,c') basis states. 
By inspection of Eq. ( 32 ) , it is easy to see that we have 



\^c'{hi,h2-Lx)). (50) 



While the states \4ic'{hi, /12)) are not single valued under 
a shift of quasihole positions by L^, path a in Fig. [6] 
can be described continuously without leaving the default 
frame. Since we have established that both \'ipc{hi, /12)) / 
and |'0c' (/ii, /12)) describe states with quasiholes in the 
same position for hi fixed and /i2 along the path a in 
Fig. |6j a relation of the form 

\iPc{huh2)) f = Y,<Ahi,h2MAhuh2)) ■ (51) 



must again hold for (some neighborhood of) this path. 
It is clear that the coefficient functions w^/(ft-i,/i2) ap- 
pearing in there must be the analytic continuation (for 
h2x > Lx) of those already defined, since 1) the argu- 



ments leading to the functional dependence Eq. ( 35 ) can 



be extended to the regime h2x > Lx and 2) for h2x < Lx 
the functions in Eq. (51 1 must be identical to those in 
Eq. (34). At the same time, for ^122: > Lx we have by 



definition 



\tpc{h2 - Lx,hi)) = '^u^^,{hi,h2)\ipc'{hi,h2 - Lx)) ■ 

^ ^ (52) 

After plugging Eqs. (|49| and (50) into Eq. (51), and 



further Eqs. (35) and (48) into both Eqs. (51) and (52), 




o- 



K: ^4 hi ^4 h 



FIG. 7. Exchange path for two quasiholes. First, the quasi- 
hole at /i2 is dragged along path Ci to ha . There the coherent 
state representation is c hanged from the original basis to the 
dual basis using Eq. ( 34 1 . The quasihole at ha is then dragged 



along C2 to hb, and the representation is changed back to the 
original basis. The quasihole at ht is moved along C3 to he- 
At this point both quasiholes are moved to their final posi- 
tions: the quasihole at h^ goes to /i2 and the quasihole at he 
goes to hi. 



comparing coefficients leads to the following additional 
relation between the f parameters: 



(53) 



All ^ parameters are thus defined up to some overall 
phase ^. We have 



V2\ 6*'^^+'"^ -1 

H- ^ A-- 



(54) 



We note that processes similar to our moves along 
global paths play a fundamental role in all studies of any- 
onic statistics on the torus (see, e.g., Ref. ,53)) . Unlike in 
the present case, it is usually assumed from the beginning 
that these anyons are entities carrying a representation 
of the braid group. Typically, complete monodromies are 
considered, where the particle moves back into its orig- 
inal position after following a path associated with one 
of the generators of the fundamental group of the torus. 
In the present case, it is of some importance that these 
global moves end before the quasihole crosses over back 
into a configuration labeled by the initial a value, thus 
changing the value of a. 



G. Braiding 



With the transition ftmctions Eq. (34) now fully de- 



fined via Eqs. ( 35 ) and ( 54 ) , the result of adiabatic trans- 
port along an exchange path as shown in Fig. [7] can be 
calculated without difficulty. We assume that in the be- 
ginning, the quasiholes are arranged at positions hi and 
h2 as shown, with h2x — hix ^ 1. The quasihole ini- 
tially at /12 is then dragged into the position he directly 
opposite the other quasihole, via path segments Ci, C2, 
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Cs which are separated by points ha, hi,. Finally, the 
quasihole at hi is moved into position /12, and the other 
quasihole is moved from /i^ into hi, completing the ex- 
change. When the one quasihole reaches the point ha, 
we pass from the coherent state expression (27) to the 



dual expression (32 1 via the transition functions, and use 



the dual coherent state expression to calculate the adia- 
batic transport along the path segment C2. At the point 
hb, the state is again re-expressed in terms of the origi- 
nal coherent state expression (27), which may be used to 



describe the completion of the exchange. 

Let the initial state be \tpc{hi,h2)), the state that lies 
in the topological sector c as defined by the Ly ^ 
limit. Adiabatic transport along the path Ci will change 
the coherent state according to 



where, using Eq. ([36|, 



(55) 



-fi=i dh'^- (i/;,(/ii, \i^c{hi, h'^)) 

JCi 



dh ■ 



(56) 



At ha we reexpress the state in the dual basis, using Eqs. 
([34| and (ISSl): 



e*^^ \Mhi,ha)) = e''''u{hi,ha)Y,Cc' \^c'{hl,ha)) . 

(57) 

We proceed by moving the same quasihole along the path 
segment C2. This process is easily described in terms of 
the dual basis states |'0c'(^i, ^2))i which appear on the 
right-hand side of Eq. ( [57| . In this basis the adiabatic 
process is simply described b y th e acquisition of a phase 
e''''^, where, using again Eq. (36), 



72 =W dh'2-{i:Ahi,h',)\V,,,JA'{hi,h',)} 
JC2 



= dh- 

JC2 



{hy,0) 



(58) 



which does not depend on the "dual" sector label c'. At 
the endpoint hi, of C2 we have thus transitioned into the 
state 



e^f^+^f^u{hi,ha)Y,Cc'HAhi,h,)) 



(59) 



The state ( 59 1 can thus be rewritten as 



e'''/^e'^'+'^'u{hi,ha)u{hh, hi)-^ <^b, hi)Cc' \i^c'{hi, ht)} 

c' 

= e^l^e'^'+'^'u{hi,ha)u{h,, /ii)- Vc(/i6, hi)) 

(61) 

The rest of the exchange path is trivially described using 
the coherent states |'0c(^i7 ^2))- The phase 73 associated 
with the path segment C3 is again given by an integral 
over a Berry connection of the form Eq. (56). The final 
move along the "baseline" C4 is carried out by moving 
both quasiholes, one from h^ into hi, and the other from 
hi into /i2. The components of the Berry connection asso- 
ciated with each complex coordinate are, however, both 



of the same form, Eq. ( 36 1 . For the remaining phases we 
thus get 



73,4 



dh ■ 



C3+C4 



(62) 



The entire exchange process thus results in the following 
transformation of the state: 



\Mhi,h2)) 



e-/2e^E.=i l^u{hl, ha)u{hb, hi)-^ \^,{hi, h2)) 

(63) 



As apparent from Eq. ( 35 ) , the u factors in the above 
equation equal i(haxhay — hi,xhi,y)/2 = —^J^^dh ■ 
(hy,hx)- When combined with the expression for 72, 
all contour integrals can be combined into a single in- 
tegral equal to the Aharonov-Bohm phase ^ab ~ Jc'^^' 
[—^{0,hx)], corresponding to a charge —1/2 particle 
moving in a unit magnetic field. We thus recover the 
well-known result [5j that the exchange of two Laughlin 
quasiparticles results in the acquisition of a phase, which 
is equal to the sum of the Aharonov-Bohm phase and a 
purely topological, statistical part: 



\Mhi,h2)) 



5 |V'c(/ll,/l2)) 



(64) 



We emphasize once more that we did not assume a pri- 
ori that any aspect of this phase is topological. Rather, 
this result followed naturally from the coherent state 
ansatz Eqs. (27), (32), and the constraints we have de- 



rived. Note that one can read the statistical phase of 
tt/2 directly off Eq. (60), which relates the transition 
functions for different quasihole configurations. While 
we have focused on the simplest case of = 1/2 for clar- 
ity, the case v = 1 /m can be treated by the same method 
through straightforward generalization [25]'|63j. 



The key observation is that this state is still in the topo- 
logical sector c as defined in the original coherent state 
basis, i.e., is of the form \tpc{hi,,hi)) times a phase. To 
see this, note that the quasiholes are now in the cr = — 1 
configuration, and we have from Eq. ( 54 ) 



^cc' ^ ^ ?cc' • 



(60) 



III. THE MOORE-READ STATE 

A. Generalized coherent state ansatz 

An appealing aspect of the method developed above, 
thus far for Laughlin states, is that the Berry connec- 
tions Eq. (36) are trivial, i.e., essentially contributing 
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only to the AB-phase. In contrast, all aspects relating 
to the statistics are manifest in the transition functions 
(cf. Eq. (60)), which need to be evaluated only at two 
isolated points. This fact might suggest that the same 
method may be amenable to discuss non-Abelian states 
in relatively simple terms as well, if suitably generalized. 
That this is so has been shown in Ref. US] for the spe- 
cial case of the Moore-Read (Pfaffian) state. In the fol- 
lowing, we will review this method, emphasizing aspects 
that need nontrivial generalization when compared to the 
Laughlin case. We will later show that the same method 
may then, with little or no further modification, be ap- 
plied to more complicated non-Abelian states also. 

The I' ~ 1 (bosonic) Moore-Read, in planar geometry, 
is the state described by the following wave function: 



■ippfizi, . . .,zn) = Pfaff 



1 



Y[iz^ - Zj) 



i<3 



(65) 

The torus degeneracy of this state is 3, and torus wave 
functions for the three ground states have been worked 
out in Ref. [M) A program similar to the one described for 
Laughlin states can now be implemented. A study 
of the special Hamiltonian [6? associated with the Pfaf- 
fian state has demonstrated that again, the three ground 
states are adiabatically connected to a thin torus limit, 
in which the ground-state patterns 111111..., 020202..., 
and 202020... emerge. 

The elementary quasihole-type excitations, which are 
again zero modes of the special Hamiltonian, turn out 
to evolve into charge 1/2 domain walls between 1111 .. . 
and 2020 . . . ground-state patterns. Periodic boundary 
conditions on the torus then require such domain walls 
to occur in even numbers. This observation is the thin 
torus statement of the well-known fact that the elemen- 
tary Pfaffian quasiholes may only be created in pairs [3]. 
For the minimum number of two quasiholes, one thus has 
four topological sectors corresponding to the sequences of 
thin torus ground-state patterns shown in Table |lj 

As in the Laughlin case, we denote these two-domain- 
wall states |ai,a2,c), and their adiabatically continued 
counterparts by \ai,a2,c). We assume that a coherent 



requirements to do this. We will thus initially assume 
(/>(/i, x) to be of the following generic form: 



(j){h, x) = exp \i(i[hy + 5/ n)x — ^{h^ — x)^ 



(67) 



Unlike in the case of the Laughlin state, we cannot derive 
Eq. ([67]) analytically from the Pfaffian 2-hole wave func- 
tions [U HI]. We observe, however, that these wave func- 
tions are holomorphic in the quasihole positions /ii, hi. 



We thus require the same for the coherent state (66), ex- 



cept for an overall normalization factor that depends on 
the quasihole positions only (and in particular does not 
depend on the parameters ai, ai in Eq. (66)). Equation 
( 67 ) is certainly the simplest expression that satisfies all 



these requirements, and is consistent with the fact that 
X and y are conjugate variables, the latter implying that 
y position enters as x momentum. The discussion of Sec. 
|II C| then makes it natural to expect that, as a function 
of /i, Eq. (67) should have the form of a LLL orbital for 



a charge 1/2 degree of freedom in a unit magnetic field 
(for some choice of vector potential, and where boundary 
conditions in hy may be twisted). This implies fi = 1/2, 
7 = 1/4, as for the i/ = 1/2 Laughlin state. We will show 
shortly that /3 = 1/2 also follows more rigorously from 
duality requirements. The parameter 7 merely controls 
the shape of the quasiholes. Its precise value will not be 
needed in the following. 

Naively, it appears that the parameter b can be for- 
mally absorbed into a shift of the coordinate origin. This 
is, however, not quite right. We will again require that 
there is a formally equivalent way to write two-hole states 
in the dual basis, defined as before via adiabatic evolution 
of domain-wall states: 



\ijc{hi,hj)) ^ Af' ^ (l){hi,Kai)(j>{hj,Ka2)\ai,a2,c) 



ai<a2 



where 



(68) 



= (const) exp [-ip{hx + S/Ti)y - j{hy - yY 



(69) 



ansatz of a form similar to Eq. ( 27 ) and its dual version It is clear that the formal equivalence between Eq. ( |66[ ) 



Eq. (32) also describe localized quasiholes in this non- 



Abelian state. In particular, we assume a Gaussian form 
for the coherent state form factors 4>{h, x) in the expres- 
sion 

\tpc{hi,h2)) = Af ^ (j){hi,Kai)(j){h2,Ka2)\ai,a2,c) 

(66) 

for quasiholes well separated along the x axis. A Gaus- 
sian form for (f>{h, x) is essentially dictated by the fact 
that X and y are conjugate variables, as argued in Sec. 



and Eq. ( 68 ) does not survive arbitrary shifts of the origin 
for the quasihole coordinates hi, /12. It is also clear that 



the coherent state expressions (66)-(69) assume definite 



II C Unlike in the case of Laughlin quasiholes, however, 
we cannot extract all the parameters entering this expres- 
sion from the analytic wave functions. Instead, we will 
have to rely more on symmetries and other consistency 



relations between the orbital indices in the LLL bases 
ipn and (fin, respectively, which define the properties of 
these orbitals under magnetic translations, Eq. ([6]), and 
determine the positions of these orbitals in space. |65| 
The choice of coordinate system, and its relation to the 
orbital indices, is also encoded in the definition of the 
symmetry operators /, r, f , Eqs. ([7|-([9|, together with 
their geometric interpretation given above. We may use 
this to severely constrain the possible values of 6. Indeed, 
these symmetries fix (5 to be a multiple of tt. Since the 
same conclusion will also emerge from duality arguments 
below, we will not pause here to show this in detail [66] . 
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The final result for the braid matrix will depend on 5 
only via e^"', which is fully determined and equals unity. 

There is one more parameter entering the generalized 
coherent state ansatz that is not yet explicit in Eqs. 
(66) and (68). This parameter enters when generalizing 



Eq. (28), which fixes the relation between the domain- 
wall positions ai.2 entering the coherent states and an 
adjacent LLL orbital with index, e.g., 2ni.2. In the case 
of Laughlin states, a single domain wall has inversion 
symmetry, and this symmetry clearly demands that the 
position a of this domain wall is defined as shown in 
Eq. (13), i.e., as the position halfway in between the ad- 



jacent ground-state patterns. More precisely, it must be 
the distance — a between this domain-wall position 
and the x position of a quasihole that suppresses the am- 
plitude in the coherent states (21) or (27). There is no 



similar symmetry argument for the Pfaffian domain-wall 
patterns. Here, quasiholes must always come in pairs, 
as mentioned above. Consider a 2-hole coherent state, 
Eq. (66), in the topological sector c = 1, Table |l] It is 



clear that the domain-wall position ai entering the co- 
herent state must be of the form oi = 2ni — s, where 
2ni is the position of the first of the string, and s 
is a shift parameter that defines the position of the do- 
main wall relative to this leading 0. For suitably chosen 
quasihole positions, an inversion symmetry leaving the 
coherent state invariant will map one quasihole onto the 
other. This does not fix the parameter s, but merely im- 
plies that the second domain wall must be assigned the 
position 02 = 2n2 -|- s, where 2^2 is the position of the 
last 0. In the topological sector c, we can thus write 



= 2ni + fi{c) , 



(70) 



where /i(l) = — s, /2(1) = +s as discussed above, and 
the values for /i(c) for c > 1 can be related to those for 
c = 1 by magnetic translations in x as shown in Table 
[Ij Here, we have defined rj = for even particle number 
iV, ?7 = 1 for TV odd. Note that the even- or oddness of 
the particle number N is just determined by the length 
of the 1111 .. . string in the patterns of Table |l] 

Equations (66)- (69), together with the shifts in the 
domain-wall positions given by Eq. ([70| and Table |l} de- 
fine the generalized coherent state ansatz. We will now 
show that this ansatz can be used to make precise state- 
ments about the statistics of the Pfaffian, and other non- 
Abelian states. 



B. The transition matrix: Constraints from 
translational symmetry 

With the generalized coherent state ansatz in place, we 
continue by carrying out steps similar to those described 
in Sees. HE and IIP for Laughlin states. Equation (34), 



the general relation between the coherent state in the 
two mutually dual bases, can be carried over unchanged. 
Again, the transition matrices appearing in these rela- 
tions are strongly constrained by translational symme- 



c 


Thin torus pattern 




He) 


1 


11111110202020201111111 


— s 


s 


2 


11111111020202020111111 


-s + 1 


s + 1 


3 


02020201 1 11111 102020202 


s - 1 


— S + T] 


4 


202020201 11111 1 10202020 


s 


-s + ri + 1 



TABLE I. Thin torus patterns for a two-domain-wall Moore- 
Read state, and the offset functions of those domain walls. 
The latter are defined in terms of the shift parameter s, and 
relate domain-wall positions ai to orbital positions 2ni (un- 
derlined) via ai — 2ni + fi{c). Some offset functions depend 
on the particle number parity rj, with -q = {rj = 1) when 
is even (odd). 



c 


T(c) 


F(c) 


1 


2 


3-^7? 


2 


1 


4: — rj 


3 


4 


2 


4 


3 


1 



TABLE II. Transformation properties of the states shown in 
Table |T] Here, it is assumed that the sector c refers to the 
original zero-mode basis, defined through the Lj, — >■ limit. 
Translating the state with would transition the state into 
sector T(c). After dragging a quasihole along the path a in 
Fig. [6] the state would transition from sector c into sector 
F{c), which is dependent on the particle number parity r;. 



try. To utilize this, we first state some of the analogues 
of Eqs. (|42|, (Esl): 



(71a) 



(71b) 

These properties again follow straightforwardly from 
the associated transformation properties of the dressed 
domain- wall states, Eqs. (39) and (40). However, the 
relation of the shifted sector T(c) to the original sector c 
is different in the present case. These relations can easily 
be read off the patterns in Table |l] and are summarized 
in Table |Tl] The remaining two transformation laws de- 
pend more critically on the value of ^, and allow us to 
determine its value. We focus on the action of Ty on 
\'4'c(hi,h2)) first. Since by duality, |-i/)c(ft.i, ^2)) is a su- 
perposition of the states |V'c'(^ij ^2)), Eq- (71b) implies 
that 



\^-)c{hi + IK, /i2 + in)) X phase factor . 

(72) 

Here, we have also used that Ty does not change the 



Ty |V'c(^l, ^2)) 



topological sector c when acting on |ai,a2,c), Eq. ( |40| . 
The left-hand side of the last equation is easily evaluated 
using Eq. (40) inside the coherent state expression. For 



1 domain-wall states, e.g., one finds 



\{ai + 0.2) for the sum in Eq. (40). With this one finds 
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that Eq. ( 72 ) indeed holds, provided that 
P = l/2, 



(73) 



Since the eigenvalues of ^t(c'),c" are doubly degenerate, 
Eq. (75) does not completely determine the coefficients 
. To this end, we must also consider the equation 



as anticipated earlier in the preceding section. With this, 
one then finds 



and similarly 



that is obtained by acting with on Eq. (34). In an 



analogous manner, this gives rise to the equation 




X!^T(c),c"^c 



(78) 



^fl fore =1,2 
1-1 for c = 3,4 



(74b) 



The relations worked out above impose strong constraints 
on the transition matrices Mcc'(^i, ^-2) defined in Eq. (34). 
We apply Ty to Eq. ([34| using Eqs. ([71b]) (with ^ = 1/2) 
and (74a). On the resulting equation, we use Eq. (34) 



again, obta ining a relation between the coherent states 

|V'c(^l, ^2)): 

X{c)e^'"^ ^ u^^, {hi + in, h2 + iK)\ipc' {hi + tK, /12 + ««)) = 



U^,,{hi,h2)e'^^''^'+'''''^/^+''\^T(c'){hl + IK, h2 + IK)) 

(75) 

where x(c) = 1 (x(c) = -1) for c = 1,2 (c = 3,4). For 
the local dependence of functions Ucc'{hi,h2) on coor- 
dinates, Eq. (35) can a gain be derived, using the same 
assuming again \hix_ 



method as in Sec. 



HE 



^2^,| 



l2y 



»1, 



3> 1. When plugged into Eq. (75), the de- 



pendence on quasihole coordinates drops out, except for 
the dependence on the quasihole configurations shown in 
Fig. [5] which is again denoted by ct = ±1. This gives the 
following equation for the coefficients , Eq. (35), 



X(c)e' 



inN—iS /-a 



T(c').c"Ccc' 



(76) 



which differs from Eq. ( 76 ) only by a replacement of the 



^ matrix by its transpose. 



To explicitly solve the constraints (76), (78), the fol- 
lowing transformation is useful. We define new topo- 
logical sector labels {fiiy), fJ,,iy ^ ±1 via the following 
superposition of states carrying c labels: 



V2 



c—3—fi/ 



(79) 



where the dependence on hi and /12 has been suppressed. 
The significance of the states lip^u^) is that under trans- 
lations in both Tx and Ty, they are now diagonal in the 
fiv label. Transition matrices u^v^^'v' and coefRcients 



^uu u'v' defined analogous to Eqs. (34) and (35), 



and are related to the quantities Ucd and via the 
transformation Eq. ( 79 ) . In terms of the matrices ^ , 



the constraints (76), (78) read 



I D = D'| D' 



(80) 



where D and D' are the diagonal matrices D = 
diag(l, 1, — 1, — 1) and D' = diag(l, — 1, 1, — 1), respec- 
tively. It is clear from Eq. (80) that only the diagonal 
elements of ^ are unconstrained, whereas the remaining 
ones must vanish. The transition matrix is thus diagonal 
in the [iv basis. We write: 



So 



(81) 



where the linear independence of the kets in Eq. ( 75 ) was 



used. For fixed c, a, this can be looked at as an eigenvalue 
problem for the quantities , c' = 1 ... 4. Obviously, 
solutions only exist if ±e*'^^~"' is an eigenvalue of the 
matrix 5t(c'),c" on the right-hand side. This is only the 
case for 



exp(2i(5) 



1 



(77) 



The coherent states are invariant, up to an unimportant 
phase, under 5 — ?> 5+211. Hence Eq. (77) narrows possible 
values of 5 down to two inequivalent possibilities. Our 
result for the statistics, however, will be the same for 
(5 = and 5 = it. We will thus keep (5 as a parameter. 



but use Eq. (77) wherever convenient. 



\i'nu{hi,h2)) = u^^{hi,h2)\%l^^,„{hi, /12)) 



£Z^u{hi,h2)\iJ^,„{hi,h2)) , 



(82) 



where u{hi, /12) is as defined below Eq. (35), and no sum- 
mation over indices is implied. We drop the tilde from 
now on, since there should be no confusion between the 
coefficient above and the coefficient defined ear- 
lier. (Note again that jj.i' should be viewed as the single 
index of a diagonal matrix element). By unitarity of the 
transition matrixes, the ^^j^'s are pure phases. 

The subscript fiv carries direct information about the 
properties of the states \ipfii^{hi, h2)) , \4'fj.iy{hi,h2)) under 
translation. From the definitions (79), it is easily verified 
directly that 
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(83) 



r 



Since Tx^ Ty are unitary operators, an expectation 
value of almo st unit mod ulus implies that the states 
\i'pLu{hi,h2)), \ipiii>{hi,h2)) are, to good approximation, 
eigenstates of these operators, with the approximate 
eigenvalue given by the expectation value. Even though 
the ^2)), \'iptj.u{hi,h2)) describe states of local- 

ized quasiholes, this is possible since and Ty trans- 
late by distances k and k, respectively, which are small 
compared to the size of the quasiholes (on the order of 
a magnetic length). To the extent that we can regard 
these states as T^, Ty eigenstates, the different associ- 
ated eigenvalues already imply that the transition func- 
tions must be diagonal in the iiiy basis, Eq. (82 1. This 
argument has been used in Ref. [25l Naively, h owever, in 
treating the states IV'/xi^C^i' ^2)), IV'm'^C^I' ^2)) as T^., Ty 
eigenstates one neglects terms that scale as 1/y/L. The 
present treatment shows that no such approximation is 
necessary in deriving Eq. ([82|). 



The transition matrix: Constraints from global 
paths 



The transition functions are thus far described by eight 
unknown phase parameters ^J^^, Eq. (82). Each of these 
parameters describes the relation between the pair of co- 
herent states {ip^i^ihi, h2)) and \tpfiiyihi,h2)) within var- 
ious patches of the two-hole configuration space. As al- 
ready discussed in Sec. 



|IIF[ these patches may be con- 
nected through paths where one quasihole is dragged 
across a frame boundary. Fig. [6j This then leads to rela- 
tions between the ^ parameters on different patches. In 
the case of the Laughlin state, all patches have been so 



connected, and there was only one independent parame- 
ter. It turns out that in the present case, the configura- 
tion space comes in two disjoint segments, which cannot 
be linked through paths as shown in Fig. [6j or any paths 
that maintain the conditions that the two quasiholes re- 
main well separated in both x and y. 

Equation ([49]) is straightforwardly generalized to the 
present case, following the same reasoning: 



\i^c{hi,h2)) f = 



ih-2yLj,+i 



''"^\^F{c){h2~LxM)) 



(84) 

Here again, / denotes a frame that will allow us to extend 
h2x beyond L^, which has been assumed in the above 



equation. Equation (49) is just a special case of Eq. 
for L = 2N + 2, 5 = TT, as befits the 1^ — 1/2 Laughlin 
2-hole state. For the = 1 Moore-Read state, however, 
one has iy = iV -I- 1 in the presence of two quasiholes. 
Also, the function F{c) assigns to c the new sector that 
one enters when the second quasihole is dragged across 
the frame boundary along the path shown in Fig. [6j The 
value of F{c) can easily be read off the patterns that 
define the 2-hole sectors in Table |l) Note however, that 
the patterns shown in the table correspond to the case of 
even particle number N, as the 1111 strings are even in 
length. As a new feature, F{c) depends on the particle 
number parity as shown in Table [Til 

Likewise, Eq. (50) may be generalized to 



|V'c'(^l,^2)) 



|Vv(^i,^2 - Lx)) 



(85) 



When analyzed in the fiv basis, Eq. (79), both the above 



equations imply that the sector labeled iiiy transitions 
into the sector labeled /i, —v when the quasihole with 
coordinate /12 is dragged along path a shown in Fig. [6] 
Specifically, Eq. (84) implies 



\i^^luihl,h2))f = 62 



^ih2yLx-\-iL5 /2 



\ipf^,-,y{h2 - Lx,hi)) X 



1 for N even, 1^ = 1, 

^^iS+ittN otherwise, 



(86) 



r 



while Eq. ( 85 ) gives 



two quasiholes in the ct = + configuration: 



\'^t,u{hi,h2)) f = i'l^^u{hi,h2)\i^ ^,u{hi,h2)) 



\^^Ahi,h2)) ^ e-"^''^^->'^\ij^,^Ah,,h2 - Lx)) . (87) 



Using the same arguments given below Eq. (51 ), we may 



apply Eq. (82) to an /-frame state \iJfj.u{hi, h2)) f with 



Here again h2x > L^, such that (hi, /12) can be taken to 
be the final configuration of the path a shown in Fig. [6] 
Plugging in Eqs. (|86| and (87) gives a relation between 



the states \i!^-^{h2 - ix, hi)) and \i}^-v{hi,h2 - L^)). 
On the other hand, these equations are, by definition, 
related via 
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\i>ti-,^ih2 - Lx,hi))f = ^f_,,^^u{h2 - Lx, hi)\%l}^,_i,{hi,h2 - L^)) . (89) 
Comparing these two relations, recalling u{hi,h2) — ei'^'^^'''^^v^^'^''^^^y\ gives the following relation between and 

f- = f+ „-^L5/2-^^M2-^.) ^ / 1 for iV even, = 1, , - 

^^'-^ ^^-^ \ ^e-«+"^ otherwise, ^ ^ 

We may also link patches of configuration space labeled by different fi, v, and a by dragging one of the quasiholes along 
path h shown in Fig. [6] This is obviously a dual version of the process just considered, and by following completely 
analogous reasoning, we find the following relation complementing Eq. (90): 



^ f+ -iLS/2-i^h{2-,.) ^ I 1 for N even, ^ = 1, , . 

f^' I jye otherwise. 



The above two equations allow us to relate any of the 
parameters with the same value of afiiy = ±1. The 
transition functions have thus been reduced to two un- 
known phases, where only the relative phase will be of 
interest. Together with the shift parameter s, this phase 
will be determined in the final step by using the locality 
considerations of Sec. IIIDI 



D. PfafHan braiding 

Given that the transition functions are diagonal in the 
fiiy basis (Eq. (82 1), the result of adiabatic exchange of 



the two quasiholes in the state \ip^}j{hi, h2)) is necessar- 
ily diagonal in this basis as well. Even in a non-Abelian 
state, it is of course possible to diagonalize any given gen- 
erator of the braid group, which describes the (counter- 
clockwise) exchange of any two quasiholes. The phase 
picked up during the exchange will, however, depend on 
the index iiv. Given the parameters defining the tran- 
sition functions, we can calculate this phase in a manner 
that is completely analogous to that discussed in Sec. 
II G In particular, the expressions ( 36 ) for the Berry con- 



nections carry over to the present case. The calculation 
is thus the same within each fii^ sector. In particular, we 
recall that the statistical part of the Berry phase could 



be directly read off Eq. ( 60 1 . Equation ( 64 1 can therefore 
be generalized to read 



^2)) 



(92) 



We denote by 7^1/ the topological part of this phase: 



(93) 



By means of the relations and (91) derived in the 



preceding section, it is clear that all phases 7^1/ can be 
related to 7++. These relations depend both on the pa- 
rameter s, as well as the particle number parity r). We 



must, therefore, distinguish the case of even (rj = 0) and 
odd (r/ = 1) particle number N. In each case, using 
i = -I- 1 we find that only even multiples of 5 enter, 
which are zero modulo 27r. Hence the parameter 5 does 
not enter the result, as anticipated earlier. For N even 
(superscript e), we find: 

7+_ = ^1+ = -7++ + 27rs, 7!L_ = 7!j_^ -f tt - Atts . 

(94) 

Likewise, for N odd (superscript o), we find: 



1+- = 1- 



-7^+ + 27rs, 7° _ = 7° + - 47rs . (95) 



There are thus three remaining parameters in the theory, 
which can be taken to be the phases 7^_|_ and 7+-(., and 
the shift parameter s. It turns out that these parameters 
are highly constrained by locality considerations of the 
kind discussed in Sec. Ill Dl 

The adiabatic transport of the quasiholes is facilitated 
by local potentials that pin the quasiholes to a certain 
location that gradually changes with time. The matrix 
elements of these local potentials in the dressed domain- 
wall basis are subject to the general considerations for 
local operators made in Sec. |IID| From these consid- 
erations it follows that the patterns contributing to the 
coherent states before and after the quasihole exchange 
can only change for orbitals whose x position (ktz, where 
n is the orbital index) is within a magnetic length (plus 
the range of the local potentials) of the exchange path. 
Regions far to the left or right of the initial quasihole 
positions do not participate in the exchange process, i.e., 
orbitals in this region are far away from any point on 
the exchange path. According to the above, this implies 
that in this region, the pattern associated with dressed 
domain-wall states entering the coherent state is unaf- 
fected during the exchange process. 

Let us consider the implications of this for the case 
where the initial state is in the sector labeled c = 3, 
Table |l] Since for a state initially in the c = 3 sector, 
all patterns form one of the two possible 2020 strings far 
to the left and far to the right of the quasiholes, this 



20 



must also be true after the exchange process, with the 
2020 patterns unchanged. This, however, imphes that 
the state is still in the c — 3 sector after the exchange. 
Identical observations can be made for the c = 4 sector. 

It is easy to translate these statements into the fiv 
basis. In order for the exchange process to be diagonal 
in the sectors c = 3 and c = 4, the phases 7^,^ must be 
independent of /i when v — ~l. This is true for both 
even and odd particle number. We thus have: 



(96) 



Note that in the case of even or odd particle number, the 
1111 strings linking domain walls in the sectors c = 3, 
c = 4 are even/odd in length, respectively. The locality 
assumptions made in Sec. II D further imply that the 



matrix elements of local operators cannot depend on the 
length of the 1111 string as long as the domain walls are 
well separated, since in this case such matrix elements do 
not depend on the separation of the quasiholes. In par- 
ticular, this implies that the Berry connection is insensi- 
tive to particle number parity (which is solely encoded in 
the length of 1111 strings) for well separated quasiholes. 
This is manifest in equations ( 36 ) which hold indepen- 
dent of particle number. However, this reasoning breaks 
down for dressed domain- wall states whose domain walls 
are not well separated. Referring to the original basis 
|oi, 02, c), this happens when two quasiholes are not well 
separated in x. In this regime, it is reasonable to ex- 
pect that matrix elements between dressed domain-wall 
states do depend on whether the (short) 1111 strings of 
patterns entering the coherent states are even or odd in 
length. This is not manifest in our formulation, since 
in this regime, we always work with the dual |ai,a2,c) 
basis. However, the transition functions that we calcu- 
lated can be expected to "know" of these parity effects. 
Hence, we expect that the phases in Eq. (96), which de- 
scribe braiding in the c = 3, 4 sectors, will depend on 
particle number parity. 

The situation is quite the opposite for the sectors c = 1 
and c — 2. Here, locality requires that the string pattern 
to the far left and right of the dressed domain- wall states 
forming the coherent states remain of the 1111 form be- 
fore and after the exchange. This only forbids transitions 
from the sectors c = 1, 2 into the sectors c = 3,4. This 
we already know from the fact that exchange processes 
are diagonal in the fiv basis, which followed from prop- 
erties under translation. However, this does not forbid 
transitions between the sectors c — 1 and c = 2. 

On the other hand, the 2020 strings forming the links 
between domain walls in these sectors, and which become 
short during the exchange process, carry no information 
about the particle number parity. This information re- 
mains hidden in the 1111 strings, which remain arbitrar- 
ily long during the exchange, in the limit of large L. We 
thus conclude that within the c = 1, 2 subspace, the braid 
matrix describing the result of the adiabatic exchange of 
the quasiholes is independent of particle number parity. 



In the fiiy basis, this leads to the following requirements: 

^ilU = g^7°+ ^ g»7l+ = f,^JU . (97) 



Using Eqs. (94) and (95), the latter reduce to the same 
equation, 7^_,_ +7°+ — ins mod 27r. Equations (96) 
give two more, 27^j_^ — 6tts—tt mod 2tt, and 27^^ 



Gtts 



mod 27r. The solutions to these equations are of the form 



3 r 

'=8-1 



= -TT TTr 

8 4 

1 1 

= -n+-7Tr, 



(98a) 
(98b) 
(98c) 



where r € Z. This amounts to eight inequivalent pos- 
sible solutions for the statistics. To discuss the relation 
between these different solutions, we first generalize our 
result to the case of 2n quasiholes on the torus. This will 
show that up to unitary transformations (taking on the 
form of simple phase conventions), all solutions are re- 
lated by overall Abelian phases. We will further obtain a 
useful pictorial representation of Pfaffian statistics, and 
relate it to more standard ones. 



E. Representation of the braid group of 2n particles 

The locality arguments used above immediately allow 
one to generalize the results obtained thus far for two 
quasiholes to the general case of 2n quasiholes. Con- 
sider the result of exchanging two quasiholes in a topo- 
logical sector as defined by taking in the — >■ limit, 
e.g. Fig. |8] Such states are the analogue of the states 
\ipc{hi, h2)) defined above, generalized to 2n quasiholes. 
Locality then implies that the result of exchanging two 
quasiholes can at most affect the string linking the as- 
sociated domain walls in the sector label. Furthermore, 
the presence of other quasiholes, which are assumed to be 
well away along the x axis, cannot affect the result of the 
exchange. One can therefore infer the result of exchang- 
ing any two quasiholes in a state of 2n quasiholes from 
the results established above for states of two quasiholes. 

These results can be generally stated as follows: 

• If the two quasiholes to be exchanged are linked 
by a 1111 string in the topological sector label, the 
state merely picks up a phase as a result of the 
exchange. This phase is given by Eq. (98b) when 



the linking 1111 string is odd in length (Fig. |8p), 
and by Eq. (98c) when the linking 1111 string is 
even in length. 

• If the two quasiholes are linked by a 2020 string, 
then upon exchange, the state will remain in the 
same topological sector with an amplitude e*^/-\/2, 
where 6 = 7r(l/8 -I- r/4 + (-1)74). It will transi- 
tion with an amplitude (— l)''ie*^/-\/2 into the sec- 
tor with the linking 2020 string shifted. [S7] 
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FIG. 8. Graphical representation of the result of exchang- 
ing two Pfaffian quasiholes for two representative pairs, a) 
A possible state in which four quasiholes could be prepared, 
labeled by its associated thin torus pattern. The state shown 
could be a four-quasihole state, in which the 20 strings at 
either end would continue around the torus, or could be a 
2n-quasihole state for n > 2, in which the ellipses mask ad- 
ditional domain walls in the thin torus pattern. The results 
of braiding any pair of quasiholes will be the same in either 
case. In the following we show only the section of the pat- 
tern relevant to the exchange; locality implies that only the 
segment of the pattern within a magnetic length of the ex- 
change path may be affected by the exchange and the rest 
remains fixed, b) Upon exchange of the indicated quasiholes, 
the state picks up the phase 7°, given by Eq. (98bl. Had the 



11 string separating the quasih oles b een even in length, the 
phase would have been 7"^, Eq. (98c I. In either case the thin 



torus pattern, and thus the topological sector of the state, 
remains unchanged, which is shown, c) When the two indi- 
cated quasiholes are exchanged the state remains in the same 
topological sector or transitions into a sector with the linking 
20 string shifted. The amplitudes for these two possibilities 
are shown next to the thin torus patterns for the sectors, 
where 6 = 7r(l/8 -I- r/4 -I- (— 1)''/4). r is an integer labeling 
the eight possible values for the overall Abelian phase, where 
r = reproduces the representation given by conformal block 
monodromies. 



These rules are represented graphically in Fig. [8] To 
make connection with the standard way to represent 
these statistics [H [TBI [171 [19] , we introduce a Majorana 
fermion degree of freedom rji associated with the i-ih do- 
main wall in the string patterns associated with our topo- 
logical sectors. Let the pair 772^, ?72j+i be associated with 
the left and right domain wall of a 1111 string. We then 
introduce fermion operators Cj = \{ri2j + irj2j+i)- Each 
Cj is now associated to a 1111 string. The topological 
Hilbert space can be constructed by acting with the op- 
erators c\ on the vacuum of the c, operators, where states 
have the j-th fermion occupied if the j-th 1111 string in 



the associated topological sector label is odd in length, 
and unoccupied otherwise. It is easy to check that ac- 
cording to the above rules, the exchange of the i-th and 
(i -|- l)-th quasihole is then represented by the operator 



(99) 



within this fermionic space, as expected for the Pfaffian 
state. The sign of 77^77^+1 in the above can be absorbed 
by a unitary transformation, facilitated by the operator 
rij ^2j - With this, the non- Abelian part of the statistics 
is thus determined unambiguously by the present formal- 
ism, whereas for the overall Abelian phase e*^, there are 
eight possible values. In the present case, these are all the 
values that are consistent with the SU {2)2 fusion rules 
[42l [43]. For r = one obtains the value that agrees 
|8] with the transformation properties of the conformal 
blocks from which the Pfaffian many-body wave functions 
are constructed [Jj . The approach discussed here is thus 
consistent with the CFT approach. For the Pfaffian case, 
the CFT approach was recently backed more rigorously 
through plasma analogy methods [10) . Similar results 
can also be obtained from the p + ip-wave superconduc- 
tor analogy [TsHTT] , although the present approach yields 
more information about the overall Abelian phase. 



IV. A: = 3 READ-REZAYI 

We have seen that the method developed above is suf- 
ficiently general to obtain the statistics of Abelian FQH 
states, and, with some adaptations, the non-Abelian 
Moore-Read state. Here we will show that the techniques 
developed in the preceding sections are indeed general 
enough to allow us, essentially without modification, to 
obtain the statistics of a more complicated non-Abelian 
state as well. We will demonstrate this for the fc = 3 
Read-Rezayi (RR) state [17] . 

Again, we focus on the bosonic "root" (highest filling 
factor, or M = 0) state of the fc = 3 sequence. This 
state has v — 3/2 and a torus degeneracy of 4. In 
taking the thin torus limit, the ground states are adi- 
abatically evolved into the patterns 0303 . . . , 3030 . . . , 
1212 . . . and 2121 . . . [35l [^. Elementary excitations 
evolve into charge 1/2 domain walls between the 3030 
and 2121 ground-state patterns, or into charge 1/2 do- 
main walls between 2121 and 1212 (Table IV). Periodic 
boundary conditions require that the former type of do- 
main wall must come in pairs, but allow the latter type 
to exist singly. Since we will need to study states with n 
quasiholes where n = 1, 2, or 3, we will begin with some 
considerations for general n. 



A. States with n quasiholes 

In the Moore-Read case, we introduced sector labels 
(/i, v) that encode the properties of states under transla- 
tions. The conventions used there made use of the fact 
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that at filling factor v = 1, and Ty commute. For 
the fc = 3 RR state at = 3/2, we thus have to proceed 
somewhat differently in exploiting translational proper- 
ties. 

To this end, we denote a thin torus state with n domain 
walls by |ai, . . . , a„; c, a), and the adiabatically evolved 
state by |ai, . . . , a„; c, a). We introduce two labels c, a 
to denote topological sectors, where a labels classes of 
sectors that are not related by translation (see Tables 
III IV), and c — 0,1 labels the two members of each 



class that are related by translation. The meaning of c is 
thus very much the same as in our discussion of Laughlin 
states. The utility of this labeling will become apparent 
shortly; the dependence of various quantities on the c la- 
bel will be constrained by translational symmetries, and 
c is conserved during braiding, in much the same way as 
for the Laughlin states. In contrast, the interesting non- 
Abelian behavior will be associated with the a label. 

We use the same mutually dual coherent state expres- 
sions as before (see Eqns. (27) and (66 1), 



E 



Y]_(t>a,jihj,Kaj) |ai,...,a„;c,a) 



ai<...<a„ j = l 



(100) 



\tpc,a{{h})) ^ Af' ^ Y[(l)a,jihj,K.aj)\ai,...,an;c,a) 

ai<...<a„ j = l 



(101) 



where the first is defined for n quasiholes that are well separated along the x axis, and the second for n quasiholes 
that are well separated along the y a xis. We have used {h} for the set of quasihole positions hi, . . . , For the same 
reasons that we discussed in Sec. Ill originally for the Pfaffian, we will assume the generic Gaussian form of ipajih, x) 
given in Eq. ( 67 ) : 



and 



'ajihj,x) = exp 



^i{hjy + S{a,j)/K)x - ^{hy^ - x^ 



exp 



'-Ah]x + S{a,j)/K)y - -f{hjy - yf 



(102) 



(103) 



In the above, we have already set ji = 1/2, which fol- 
lows in exactly the same way as for the Pfaffian. We 
have written ^ as a function of the sector a, to allow for 
the possibility that the momentum shift 5 may take on 
different values for quasiholes associated with different 
types of domain walls. However, (j) is independent of c 
since the type of the j-th domain wall is invariant under 
translation. 



Again, the two bases ( 100 1 and ( 101 ) are related to 



each other by a transition matrix. In general, the ele- 
ments of this matrix depend on both c and a. 

c' ,a' 



In complete analogy with Eq. (35), we can derive the 



local dependence of the transition matrix within each of 
the regions labeled by a, which are components of the 
quasihole configuration space with quasihole coordinates 
well separated in both x and y (cf. Fig. [s] as well as Fig. [9] 
below). 



c'a,.'(W)-C.c'aW«({M), 



(105) 



again with u{{K\) = ^j'^i^'^iy _ Por n — 2, there are 



72 of these parameters f^^, ^ we distinguish two con- 
figurations a (Fig. [5| , and for each there is a 6 x 6 matrix 
in the sector labels. 

We first state the translational properties of the n- 
domain-wall states, which are the same as in Eqs. (39) 
and (40 1, since a is a spectator under translations. We 
now adopt a natural definition for the c labels. Recall 
that the action of Ty is given as follows. 



Ty |ai, . . . , a„; c,a) = exp 



2m 



|ai, . . . ,a„; c, a) 
(106) 

where the Uj are the orbitals occupied in the pattern 
labeling the state. We find that the sum over the nj 
takes on the following form. 



(107) 



where c = 0, 1, and the domain- wall positions are defined 
via flj ^ 2ni + /i(c, a) as before, with the orbital position 
2ni defined in relation to the domain wall as shown in 
Table IV Equation ( |107 ) then defines c modulo 2, and 
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a labels the three "supersectors" formed by the transla- 
tional pairs of states. 

The translational properties of the n-domain-wall 
states then are 

Tx\ai,. . . ,a„;c, a) = |ai + 1, . . . , a„ + 1; 1 - c, a) 
:r |ai, . . . ,a„;c,a) = e ^ 3 j^^^ 



Ty \ai, . . . , a„; c, a) 



(108) 



|ai, . . . , a„; c, a) = |ai + 1, . . . , a„ + 1; 1 - c, a) 



(109) 



where again we write A = lyL, this time with v — 3/2. 

As in the preceding cases, the translational properties 
of the coherent states follow directly from Eqs. ( 108 ) and 



(1091 



and similarly used {/i + zk}. We can 
use these translational properties to completely deter- 
mine the dependence of the transition matrices on c. To 
make this decoupling more explicit, we introduce two- 
component objects, denoted by a 



l^a(W)) = 



l*a({M)) 



IV'0,a({M)) 



(112) 



7;|Vc,a({M)> - e^"^"''"1^c,o({/i + 4)) (110) 

r,|Vc,a({M)) = e^-E, + *("-J")|^i_,.,({/i + z«;})) , 

(111) 

where we have used the notation {/i + k} = /ii + 



where the phase splitting between the c = and c = 1 
states has been introduced with foresight to keep later 
phases in check. Correspondingly, we may view the full 
transition matrix as a "supermatrix" 2'^, i.e., an amax x 
Climax matrix, the elements of which are each 2x2 matrices 
denoted ^, . So we write the equation between the 
original and dual bases as 



i*a(W)>-E^(W)^s."'i*"'(W)) 



(113) 



Note the similarity of Eq. ( 113 1 to the Laughlin transition 
matrix Eq. (44), to which Eq. (113) reduces for amax = l- 



Wc rewrite Eqs. (110) and (|111|) in terms of the two- 
component basis. 



Ty\^^m))^^~ 



1 

-1 



|*a({/i + *«})) ; T,|*„({/i})) - e^"^^. 







1 



(114) 



|*„({/i + ik})) 



(115) 



Eq. (113), Eqs. (114) and (1151 each give a consistency 



As before, (cf. Eqs. (76), (78)), when applied to together with the constraint 

27riA -|- i (5(q:, j) 



equation that must be satisfied by every ^, : 



exp 



which imply 



g*Ej S{a,j) 



1 \ 



1 

-1 







-1 j I e^5:,^("j) 



'a.a' 



1 e"^ 



(118) 



(119) 



(116) In the above, is an overall coefficient, and y/2 is 

normalization factor. 



The phase choice we made in Eq. (112) has allowed 
us to decouple the a and c indices within the transition 
function. We can write the matrix S"^ defining the tran- 
(117) sition function Eq. (113) as 

E'^ (^M , (120) 



where is the amax x Omax matrix of coefficients ^, 
and M is the 2x2 matrix 



M = 



1 
71 



1 e-'""^ 



(121) 
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Sector c, a 



0,1 
1,1 



Thin torus pattern 



IVIVIVIIVIVIVIVIX 
VIVIVIIVIVIVIVIVI 



hie, a) 



TABLE III. Thin torus patterns for a single-quasihole k=3 
Read-Rezayi state, and the offset functions of the associated 
domain walls. The latter are fully determined by inversion 
symmetry of the state. The orbital positions, 2ni, are under- 
lined. Since the sectors are all related by translation, a takes 
on a single value. 



The a dependence of S*^ is completely contained in the 
corresponding coefficient matrix and the c depen- 
dence is completely contained in the M matrix. 

If we consider the translational properties of the states 
in the case of a single quasihole, we can constrain some 
of the S{a,j) parameters appearing above. For a single 
quasihole on a torus, the only topological sectors respect- 
ing periodic boundary conditions are those with domain 
walls between 2121 patterns, as shown in Table [TTTl There 
are two such sectors, related by translation, so for a sin- 
gle quasihole amax = 1- There is then only a single S {a, j) 
parameter, which we call d. When we consider Eq. (119) 
and note that in this case A — vL = |(|(2iV 
half-odd integral, we find d — ir. 



iTTis 



In general, the 5(a, j)s are each associated with a cer- 
tain type of domain wall, so by fixing d in the single- 
quasihole case we also fix any S{a,j) associated with a 
212l|l212-type domain wall in an n-quasihole state. We 
can constrain the other S{a,j)s to be either or tt by 
considering Eq. (119) in the case n = 2. For two quasi- 



holcs there are only two independent S{a,j) parameters: 
5(3, j), which is associated with 2121jl212-type domain 
walls and is thus known to be tt from the one-quasihole 
argument; and S{l,j) and (5(2, j), which are associated 
variously with domain walls between 1212 and 0303 
strings, and which must be equal by the argument in Ref. 
I66l For n = 2, we have A = i/L = |(|(27V-h2)), which is 
an integer, and Eq. (|119l) reduces to exp[i d{a, j)] = 1. 
This is already satished for S{3,j) — tt, and can be sat- 
isfied for a = 2,3 only if (5(1, j) = (5(2, j) = 0,7r. 

In the end, we want to find explicit expressions for the 
elements of the transition matrices S"', which we have 
reduced to the problem of finding the elements of the 
coefficient matrices. This will be our task in the following 
sections. 



B. Two quasiHoles 

The thin torus patterns for two-quasihole states with 
c = are given in Table |IV[ To find the statistics of these 
quasiholes we must constrain the transition matrices S"*" 
and S~. Both S'^s have nine complex unknowns, the en- 
tries of the ^'^ matrices. To constrain these we will move 
the quasiholes around global paths, which we defined in 
Sec. IIIFI We will then make further use of the mirror 
Thin torus pattern fi{<^) f^iot) F{a) 



303030302121212030303030 
121212120303030212121212 
121212112121212112121212 



1 - s 



1 + s 



TABLE IV. c = thin torus patterns for a two-quasihole 
fc = 3 Read-Rezayi state, and the offset functions of the as- 
sociated domain walls. The orbital positions, 2ni, are under- 
lined. Patterns for c = 1 can be obtained by shifting each 
occupancy number one orbital to the right, and c — 1 offset 
functions by adding or subtracting 1 to each offset function 
above, whichever is more convenient. 



symmetry operation, which has thus far only been dis- 
cussed in Sec. |I A| and very briefly in Sec. |III B[ As in 
the Moore-Read case, we gain further constraints by im- 
posing locality and unitarity. In the general solution to 
these equations some unknown parameters still remain. 
We will be able to constrain the latter by subsequently 
studying the case of three quasiholes in Sec. |IVC| 



Constraints from global paths 



As discussed in sections II F and IIIB the transition 



matrices for different configurations can be connected 
by dragging the quasiholes through the global paths in 
Fig. |6] We first consider two quasiholes in the cr = + 
configuration, and imagine the right quasihole moving 
around the x direction of the torus along the path a in 
Fig. [6j Using the reasoning of Sec. |III C| we find the 
following effects on the coherent states: 

(122) 

Moving the quasihole along this path changes the sector 
label a for the original basis into F{a), the values of 
which can be read off the patterns and are summarized 
in Table HVI 



To find a constraint on the ^°'s, we write Eq. ( 122 ) in 
the two-component basis. 



25 



1 

-1 



g-5iEj'5("j) 



where we have used that fj{a) = fj{0,a) = fj{c,a) — c mod 2. Applying Eq. (123) to Eq. (1131 gives 



a' \ 

We can simpUfy Eq. (|124| using Eq. ([Tl6|). 







1 

-1 



\'^F{a){h2 - L^, hi)) = u{h2 - L^, hi)e 



A-,:-|5(a,2)-i*i:,'5("j)"+ 



-j7r/2(a') I 



(123) 



(124) 



(125) 



We want to write this as an equation between and 
which we can do by noting the equivalence between 
Eq. ( 125 ) as written and Eq. \\12>\ evaluated at quasihole 
positions (h2—Lx, hi). To make this equivalence manifest 
we can write the action of F in matrix form as: 



I 

inverted to yield the following equation: 












A 




p 







A 














I 








iB)c 



B = 



'Ja,F{a'} 




(126) 



(127) 



(130) 



Combining Eqs. (129 1 and (1301 gives us a nontrivial 
consistency relation for 




Ap 

Ap-i 
0-1 



(131) 



Since the transition matrix in Eq. ( 113 ) evaluated at posi- 
tions (/i 2 — L x, hi) involves S^, and the transition matrix 
in Eq. (125) is a product involving S+, the equivalence 
of these two equations implies: 



C = B^Miag 



^+diag 



T/2(a) 



(1^28) 

where we canceled the matrix M common to both S'^s, 
and the argument of diag[. . . ] specifies the a-th diagonal 
entry of a diagonal matrix. If we use the values of /2(a) 
from Table Irvl Eq. becomes 



Equation (131) gives us several equations between the 
matrix elements of the coefficients S.^^,- Equation 
(129 1 reduces the number of unknown ^,s from eigh- 
teen to nine. The consistency relationship Eq. (131) fur- 
ther reduces the number of unknown elements from nine 
down to five. A particular choice for the five independent 
iaa'^ is the following: 



62 P^Cii 
C31 -AK31 




(132) 





^0 


A 


0^ 






A 












^0 










p 


0-1 



(129) 



We have defined two phases: p = — exp [i7r(s -I- |)] and 
A = exp [i(L/2 -I- l)(7r — (5)]. Note that for two quasi- 
holes L is even and A^ = 1. 

We can perform the same process in the y direction 
and drag the quasihole around the global path marked 
b in Fig. [6] After an argument similar to that above we 
find another equation between ^~ and which can be 



Note that any time the configuration index cr is omitted 
as in the above equation, we take it to be (t = +. 



2. Constraints from mirror symmetry 

We now make use of the antilinear symmetry opera- 
tor r defined in Sec. |I A[ i.e., the combination of time 
reversal and mirror symmetry. Applying r will exchange 
the X positions of the quasiholes across the y axis. This 
operation changes the configuration a, which will allow 
us to derive another equation between ^+ and First, 
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we describe how this symmetry acts on an ri-quasihole 
state. 



states, thus its action on those states is: 

T |ai, . . . , a„; c, a) = \L - a„, . . . , L - ai; c, F^(a)) 



From the definition Eq. ([S]), the effect of r on bare LL 
product states is clear: it reflects the original basis states 
across the y axis, and it has no effect on the dual states. 
For bare product states with domain walls, the domain- 
wall positions will be similarly reflected, r commutes 
with the adiabatic evolution operators (as constructed, 
e.g., in Ref. that define the delocalized quasihole 



r|ai. 



,a„;c, a) = \ai 



(133a) 
(133b) 



We write that the position of the j'-th dual-basis quasi- 



hole aj goes to L — aj in Eq. ( 133a ) so as to stay within 



the default frame. Also note that in general r might or 
might not change a, and we describe this change by some 
function Fr, the values of which can be found from the 
patterns. It turns out that for the case of two quasiholes, 
i<V(a) — OL- Later when we analyze the case of three 
quasiholes, will be a nontrivial mapping. 

Equation (1331 allows us to derive how r acts on co- 
herent states of n quasiholes. In terms of two-component 
states: 



T|*a({/l})) = (")|*a({-/l* +i,})) . 



(134a) 
(134b) 



For now, we will restrict ourselves to the case of two 



quasiholes. In this case, Eq. (134) simplifies to 



T \^o.{hi, h2)) = E, h,, (135a) 

(135b) 



\^c.{hl,h2))^\^c{K,h',)), 



where for all indices 



-h* 



To arrive at 



over from Sec. |II G| because all the Berry connections 
along the path segments considered above are indepen- 
dent of the sector (see Eq. (56 1 for example). I.e., for 
the exchange of two quasiholes as in Fig. [7) dragging the 
second quasihole along the path segment Ci causes the 
wave functions in each sector to pick up the same phase 
exp [171] defined in Eq. (55). 



Eq. ( 135 ) we have used that for two quasiholes the phase 
factors on Eq. (134a) and Eq. (134b]) are both 1 — the 
former because L is even, and the latter can be seen by 
inserting the values of fj{a) from Table [Tv] — and Fr{a) = 
a as noted above. Equation ( 135 ) allows us to apply r 



/ \^i{hi,h2)) 
\^2{hiM)) 
V \^i{hiM)) 



( |*l(/il,/ia)) 

\^2{huK)) 
\\^3{huK)) 



(137) 



Reiterating the remaining steps described in Sec. II G 



to Eq. ([113f . Let us begin with the two quasiholesin the ^^^^ ^^g^j^ adiabatic exchange is the following: 



a — + configuration; when we apply t to Eq. (113) and 



compare the resulting equation to Eq. (113) evaluated 



at the changed spatial coordinates, we find the simple 
relationship r, = (Ca')*> or 



|*l(/^l,/i2)) 
|*2(/il,/i2)) 



-> e 



\^l{hi,h2)) 
|*2(/il,/i2)) 
|*3(/ll,/l2)) 



(136) 



For the moment, we leave the relation ( 136 1 implicit, and 
use it in App. |X] to further reduce the number of inde- 
pendent parameters. 



3. Braiding 

We can perform the adiabatic exchange of two quasi- 
holes using again the method of Sees. |II G| and |III D| 
with minor generalizations. The details formally carry 



(138) 

Once again we see that adiabatic exchange results in a 
path-dependent Aharonov-Bohm phase and a topologi- 
cal, statistical part made of a product of the transition 
functions, which we call the braid matrix. The structure 
of this matrix is 



X' 



112x2 , 



(139) 



where the translational, c-dependent part of the braid 
matrix is the product MM'^ = I2x2, and we have defined 
the "reduced" braid matrix as the a-dependent part, 



(140) 
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We can constrain the form of the matrix x by making 
an argument about the locahty of the exchange process, 
analogous to the argument made in Sec. |IIID| Recall 
that according to the latter, only the string of the pat- 
tern that is between the domain walls taking part in the 
exchange can be changed as a result of this process. Any 
regions of the pattern far to the left or right of the ini- 
tial positions must remain unchanged after the exchange. 
For one, this requires the exchange processes to be diag- 
onal in c. This is already manifest by the structure of 
the braid matrix derived thus far, Eq. (1391. However, 



certain transitions of the a label are allowed. To see this. 



we again refer to Table IV One observes that transitions 
into and out of the a = 1 sector are forbidden, since this 
is the only sector with 3030-type patterns far to the left 
and far to the right of the domain walls. The other two 
sectors have 2121-type patterns at the left and right end. 
Therefore, transitions between these sectors are allowed. 

These considerations imply that the reduced braid ma- 
trix Eq. ( 140 1 must be of the form: 







(141) 



where dots indicate (potentially) non-zero matrix ele- 
ments. Equation ( 141 ) gives two independent constraint 
equations for the matrix elements £,a,a' ■ We will also use 
constraints for the ^q,q'S gained from the fact that ^+ 
must be unitary: 



e+(e+)^-i 



3x3 



(142) 



Together, Eq. ( 141 ) and Eq. ( 142 1 provide enough con- 
straint equations to fix the S,a.a' up to the parameter p, 
introduced after Eq. ( 129[ ), which is in turn defined by the 
shift parameter s defined in Table [TV] These constraints 
allow us to write explicit expressions for the elements of 
the braid matrix in terms of only the parameter p: 



( 



\ ±v?+jFi(i-p) p + 





P^ii-p) 

0-1-1 



(143) 

The details are presented in Appendix [A] While p is still 
unknown at this stage, it is no longer completely uncon- 
strained. To further constrain the value of p and fully 
determine the statistics, we must study the case of three 
quasiholes. 



C. Three quasiHoles 

We expect that we can gain new information about the 
statistics by braiding two quasiholes among a system of 
three. To see this, note that as long as there are only two 
quasiholes, boundary conditions require that both are as- 
sociated with the same "domain- wall type". I.e., both 



a 


Thin torus pattern 


/i 


/2 


h 


F 




1 


30303021212112121203030 


s 


1 
2 


1 - s 


3 


1 


2 


12121121212030303021212 


1 

2 


— S 


s 


1 


3 


3 


12121203030212121121212 


1 - S 


-l + s 


1 

2 


2 


2 


4 


12121121212112121121212 


1 
2 


1 
2 


1 
2 


4 


4 



TABLE V. c = thin torus patterns for a three-domain- 
wall fc = 3 Read-Rezayi state, and the offset functions of 
those domain walls. The orbital positions, 2?!^, are under- 
lined. Patterns and offset functions for c — 1 can be obtained 
by, respectively, shifting each pattern one orbital to the right 
and adding 1 to each offset function. 



domain walls must either occur between a 3030 string 
and a 2121 string, or between two 2121 strings. Hence, 
we were not yet able to study what happens when a quasi- 
hole associated with the former type is exchanged with 
one associated with the latter type. To study such pro- 
cesses, we must consider systems with three quasiholes. 
The relevant topological sectors are displayed in Table [Vl 
It will suffice to exchange the first two quasiholes (along 
x). The "new" situation described above will then occur 
in the sectors a = 1 and a = 2. Using locality arguments 
analogous to the preceding section, we conclude that ex- 
changing the first two quasiholes in these sectors is a 
diagonal process, since it is not possible to reach a differ- 
ent sector by replacing the string linking the associated 
domain walls. On the other hand, by complete analogy 
with the preceding section, the same exchange processes 
may lead to transitions between the a = 3 and a = 4 
sectors. These processes are locally the same as those 
discussed for the a = 2 and a = 3 sectors in the preced- 
ing section. Invoking again locality, within the a = 3, 4 
subspace the (reduced) braid matrix must be given by the 
same 2x2 block displayed in Eq. ( 143 ) . We used exactly 
the same argument before in Sec. 



HI E where we con- 



structed the 2n-quasiparticle representation of the braid 
group from the two-quasiparticle braid matrix for the 
Moore-Read state. These arguments constrain the form 
of the reduced braid matrix associated with the first two 
quasiholes to be: 



V 



p{p + p 
±^/p + p^{l-p) 



1) ±^p+p-^{l-p) 



p + p-^-l J 

(144) 

where the dots indicate some matrix element we do not 
yet know, and blank spaces represent zeros. In the above. 



p is the same parameter appearing in Eq. (143), but we 
leave it understood that the quantities x, ^ and S in this 
section refer to the three-quasihole case, and arc different 
from their two-quasihole counterparts. In the above, we 
have anticipated that braiding will again be diagonal in 
the c label, and x is again defined through the action of 
braiding on the a label, which will follow below. 

We may again proceed by expressing x through the 
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FIG. 9. Configurations a = {cj\,02)- (Ji indicates the relative 
position of the leftmost two quasiholes, (J2 indicates the posi- 
tion of the third quasihole relative to the first two. Top (left 

to right): +- 1-, +0, and H — . Bottom (left to right): , —0, 

and — h. 



transition matrix coefRcients ^, and deriving various 
constraints on the latter, where now additional con- 
straints come from the 2x2 block in (144). The pro- 



cedure is analogous to the preceding section, where only 
one aspect requires nontrivial generalization: in the two- 
quasihole section there were only two transition matrices, 
S"*" and S~, one for each configuration. For n quasiholes, 
we must distinguish n\ configurations and define a tran- 
sition matrix for each. We choose the following notation 
to label these configurations. For an n-quasihole system, 
we let a ~ (ui, . . . , cr„_i). cti takes a value -I- or — , 
indicating the relative position of the two leftmost quasi- 
holes, in the same manner as in the preceding section. 
(72 takes a value -f , — , or and indicates the position of 
the third quasihole relative to the first two, as shown in 
Fig. [9] for three quasiholes. We could proceed further in 
the same way for n > 3 quasiholes, but n < 3 suffices for 
our purposes. 

With these conventions, we can find the result of ex- 
change of two quasiholes in terms of the transition matri- 
ces. As pointed out, we choose to braid the two leftmost 
quasiholes and leave the third fixed. Further, let us say 
it is fixed "above" the other two, so (T2 = +. |M| The 
exchange can be broken down into segments in complete 
analogy with the t wo-q uasihole case, yielding an equation 
analogous to Eq. (138 1: 



l*2({M)> 
l*3({M)) 
V l«'4({M)) / 



*2({M)) 
*3({M)) 
V l*4({M)) J 

(145) 



As anticipated above, using Eqs. (120), (121), we find 

S++(S"+)t - 



X' 



112x2 



(146) 



which again defines the reduced braid matrix x in terms 
of the coefficient matrices, 



x = r+(r+) 



(147) 



As in previous examples, we will use symmetries and 
global paths to constrain the ^'^ matrices, then use the 
implications of locality, Eq. ( 144 ) , to find explicit expres- 
sions for the elements of x- 



1. Constraints from, mirror symmetry 



In Sec. |IVB2[ the action of T on n quasihole states 
has been discussed, Eq. (134). The three-domain- wall 
patterns are shown in Table |Vj along with the values of 
Ft- (a), which follow directly from these patterns. We can 
represent the map -FV(a) in matrix form: 



Br = 



/I 



V 



1 

1 



1/ 



(148) 



If we apply r to Eq. ( 113 ), proceeding as in the derivation 
of Eq. ( 136 ) and using information from Table [v] we find 



/ A2 



V 







A2 




in* 



1/ 



1 



1/ 



where the phase A = 
odd, so A^ — exp[i(7r 



exp[-i^(7r - 5) 



(149) 
. Note that L is 
1. The function 



5)] and A-* 

grio') gives the new configuration after reflection of the 
quasiholes in configuration a across the y axis, and its 
values are given in Table |VI| 

The matrix structure of the last equation is somewhat 
more complicated than Eq. (136). Unlike the latter. 



Eq. (149) is not "self-dual", i.e., we may obtain an anal- 
ogous but different equation by using the "dual" mirror 
symmetry operator f instead (Sec. I A). It reads 





/ 


1 




1 




V 



/ A2 



A2 
A2 



1/ 



(150) 



The function 5t(o') captures the change in configuration 
under f. Its values are given in Table [Vl} 



We now evaluate Eq. (149) for a — ( — ,—), ,9r = 
{+, +), and Eq. ([150]) for a = (+, +), g^{+, +) = (-, -), 
and plug one into the other. This gives the following 
consistency equation for ^"'"^: 
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A2 
A2 



A2 
A2 
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+0 
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+0 


-+ 


-+ 


-0 


+- 


-+ 
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-0 



TABLE VI. The effect of various operations on the three- 
quasihole configuration a. g{cr) is the resulting configuration 
when a quasihole is dragged around the torus along an a- 
type path (as shown in Fig. |6] for two quasiholes), and g(cr) 
is the same for a fe-type path. The rightmost and topmost 
quasihole is being dragged, respectively, ^^(f") is the resultant 
configuration under the mirror refiection r, ^((j) is similarly 
defined for t. 



which constrains to be of the form 
f en 



V ^41 



31 



^23 



^13 
^23 
^22 
^43 



Cl4 \ 

34 
^34 
C44 / 
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As before, when the configuration a is omitted, we take 
it to be ++. 



2. Constraints from global paths 

We continue with our program by deriving constraints 
from "global paths", as done for the two-quasihole case 



in Sec. |IVB 1 We begin by generalizing Eq. ( 123 ) 
(cf. Fig. |6| to the case of three quasiholes. In the two- 
quasihole case we assumed the two quasiholes to be in a 
a = + configuration, then moved the top right quasihole 
around the x direction of the torus to the top left. In this 
section we will need to derive more general behavior, al- 
lowing that the rightmost quasihole can be at the top, 
middle, or bottom relative to the other two quasiholes. 
The analogue of Eq. ([l22| 
in Fig. [6) then becomes: 



for a path similar to path a 



(153) 



Here fj{c,a) can be inferred from Table |vj and j equals 
1, 2, or 3 if the quasihole encircling the torus is respec- 
tively the first, second or third when viewed from the y 
direction. The position {h'} = — Lx,hi,h2- As be- 
fore, the change in a after moving the quasihole along 
the path is described by the function F{a). Its values 
directly follow from the associated patterns, as discussed 
m S ec. |IVBH and they are given m Table |V| We recast 
Eq. ( 153 ) in the two-component basis: 



l*a({M)) 



(154) 



1 



l*a(W)) = e-"-^^("' I ^ )l*a(m)) 



Just as in the two-quasihole case, Eq. ( 154 ) allows us it has the form 



to derive an equation between the transition matrix in 
the configuration a with the matrix in the configuration 



B ^diag[e 



^5(a,3)-i«E,«("j)]^-diag 



-i-irf j{a) 



(155) 

The pairs {a,g{a)) are summarized in Table VI The ma- 



trix B is defined as in Eq. ( 126 1, and for three quasiholes 



B = 



( 1 \ 

1 

1 

V 1/ 



(156) 



In any specific instance of Eq. (155), one first chooses a 



starting configuration a, and identifies the corresponding 
y-direction index j of the quasihole that will encircle the 
torus, j is in one-to-one correspondence with a2'- for a2 
is +, 0, or — , j is respectively 3, 2, or 1. For instance, 
were we to begin in configuration -|— 1-, then j = 3 and 
after the encircling the system would be in configuration 



30 



-0. Thus we find the relation 



A 

A2 



A 



impHed by locahty, as we argued above. The details are 
given in App. 151 resulting in the following form for x'- 



1/ 



-1 



/ 



-1/ 



where p is defined as before, p 



(157) 



p\l-p) 



\ e'-'^'^p^^/p + p-^ - 1 



e^^'^^p{l-p) 



p-^-l 



J 



We can go through the same derivation for a path sim- 
ilar to path b in Fig. [6j in which the top quasihole moves 
around the y direction of the torus and ends at the bot- 
tom. We find 



(162) 

where 6 and 62 are as yet undetermined phases. 

In deriving the above equation, only the zero matrix 
elements of Eq. ( 144 ) have been used. To enforce consis- 



ts ('^^ = diag 



tency between the two- and three-quasihole braiding ma- 
t'^diag[e~'''^~*^'''^"'^''~5»Ii:j trices, as dict ated by locality, we must equate the 2x2 

block of Eq. (fT62l) to that of Eq. plif. Equating the 



(158) 

It is very important to note that for this path, the mean- 
ing of the index j is different from the previous path. In 
the previous case, the quasihole encircled the torus in the 
X direction, so the moving quasihole was the rightmost in 
horizontal (x) order but was the j-th quasihole in vertical 
(jj) order; in this case, the quasihole encircles the torus in 
the y direction, so the moving quasihole is the topmost 
in vertical order but is the j-th in horizontal order. As 
an example, if we begin in configuration —0, the topmost 
quasihole is that on the left, so j = l. Plugging in the 
appropriate values from the tables. 



fp- 



V 



A 
A2 

-1/ V 1/ 



(159) 

Combining Eqs. (1591 and (157) gives us a consistency 
relation for t^"*". 



/ Ap- 
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Ap 
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Ap-i 
-A2 



which further constrains in addition to Eq. ( 152 ) 



/ ar 

V ^41 



13 



A^e 
AV61 
~A-va3 

-A-^pUi 



63 

-A-ipCi3 
-ApUi 



64 

-A-1K14 
-AK14 

^44 
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3. Braid matrix 



expressions for the element X33 gives us a consistency 
relation that we can use to constrain p: 



(163) 



If we define x=p + p ^ for convenience and take the ab- 



solute square of Eq. (163), we find 

2 - x = (a; - 1)2 , 

which is solved when x is the golden ratio, 

l + %/5 
X = ip — — - — . 



(164) 



(165) 



We have chosen the positive root because Eq. ( B9b ) im- 
plies X > 1. If we define the angle a by p — exp[iTTa\ then 
x = 2cos(7ra) and Eq. (165) implies 



(166) 



n is also related to the shift parameter s = a -1- 1/2, and 
sdEq. ([166]) tells us 69J 



-1/ 



-^4 



(167) 



n 60) The phase information in Eq. ( 163 ) fixes the overall phase 



e'^ = e"" . (168) 

There are two more consistency equations found from 
equating Eqs. (144) and ( |162[ ). One yields exp [z02] = 
±1, and the other is trivially satisfied when x ~ ip. Up 
to some signs, the braid matrices for two- and three- 
quasihole systems have thus been completely solved for. 
We will discuss our solution(s) in the following section. 



As in the two-quasihole section, we will further deter- 
mine the structure of the reduced braid matrix using con- 
straint equations from unitarity and from locality. En- 
forcing locality means that we equate the matrix product 
for X in Eq. (145) with the form in Eq. (144), which is 



V. DISCUSSION 

In Sec. |IV[ we have found solutions for the braid matri- 
ces describing exchange processes between two and three 
quasiholes that are consistent with the coherent state 



31 



and the three-quasihole matrix is 




.30302121203030302121211212120303(£ 
^ 



...303021212030303021212112121203030.. 
...303021212030303021212112121203030.. 



+ 



X e 



...303021212112121121212112121203030., 



...303021212030303021212112121203030... 



FIG. 10. Graphical representation of the resuh of exchanging 
two fc = 3 Read-Rezayi quasiholes for three example pairs. 
Top) A possible state in which five quasiholes could be pre- 
pared, labeled by its associated thin torus pattern. The state 
shown could be a five-quasihole state, in which the 30 strings 
at either end would continue around the torus, or could be an 
n-quasihole state for n > 5, in which the ellipses mask addi- 
tional domain walls in the thin torus pattern. The results of 
braiding any pair of quasiholes shown here will be the same 
in either case, a) Upon exchange of the indicated quasiholes, 
the state picks up the phase e"^'!"")^ where a is given by 



Eq. ( |166[ ) (with the lower sign correctly describing the con- 
formal block monodromies of the RR trial states) . The thin 
torus pattern, and thus the topological sector of the state, re- 
mains unchanged after the exchange, as shown, b) When the 
two indicated quasiholes are exchanged the state remains in 
the same topological sector or transitions into a sector with 
the linking 30 string changed to a 21 string. The amplitudes 
for these two possibilities are shown beneath the thin torus 
patterns for the sectors, where is the golden ratio, Eq. ( 165 1, 
c) Upon exchange of the indicated quasiholes, the state picks 
up the phase e"'^2+2''). The thin torus pattern, and thus the 
topological sector of the state, remains unchanged after the 
exchange, as shown. 



V 



I 



(170) 

In vifiiting these matrices, we have removed the ± from 
the off-diagonal elements, choosing the -I- sign. Choosing 
the negative sign instead leads to a unitarily equivalent 
solution, where the transformation is facilitated through 
multiplication of each state by (— 1)^^^°, where #30 is 
the number of 3030 . . . strings in the thin torus pattern 
associated with that state. The arguments given below 
will make it obvious that this equivalence also carries over 
to general n-quasihole sectors. We have thus obtained 
only two unitarily inequivalent solutions. It is clear from 
the above that these two solutions are closely related, 
namely by complex conjugation and an overall Abelian 
phase —1. Thus, the non- Abelian content of the fc = 3 
state has been determined uniquely by our method. 

We will now use the locality arguments already made 
in Sec. |1V B 3| for states of three quasiholes, and applied 
earlier in Sec. |I11 E| to the Pfaffian case, to generalize 
these solutions to the case of n quasiholes. In essence, 
these arguments implied that the result of exchanging 
two neighboring quasiholes can only affect the ground- 
state pattern linking the associated domain walls in the 
sector label, and only depend on the sequence of three 
patterns that are separated by these two domain walls. 
For this, however, all possibilities have been exhausted by 
considering two and three quasiholes, respectively. We 
can thus list the following rules, applicable to general n- 
quasi hole states, obtained directly from Eqs. (169) and 
( [T70| : 

• If the two quasiholes to be exchanged are associated 

with domain walls between ground-state strings 

. . . 303d212lll21 ... or . . . 212lll212l030 . . . , then af- 
II I I ' 

ter exchange the state remains in the same sector 
and picks up the phase e"^(2+2«). 



ansatz for the fc = 3, j/ = 3/2 Read-Rezayi state. In the 
following, we discuss how many independent solutions we 
have found, how they lead to general rules for the braid- 
ing of n quasiholes, and how these solutions compare to 
those obtained by other methods. 



By means of Eq. ( 165 1 , we may now express the braid 
matrices for two quasiholes, Eq. (143), and three quasi- 



holes, Eq. (162), in terms of only the golden ratio (/s and 
the parameter a = ± ^ . The two-quasihole matrix is then 



/ 



(169) 



• If the two quasiholes are associated with the pat- 
tern . . . 3030l2121^030 . . . , then after the exchange 
the state merely picks up the phase e"'-'"'^). 

• If the quasiholes are associated with the pat- 
tern . . . 212|0303d21 . . . , after exchange the state 
will stay in same topological sector with am- 
plitude e"'^2"'""-*(/?~^ s-nd transition with ampli- 
tude e"*-2~^°)(^~ 2 into a sector that has the ex- 
changed quasiholes associated with the pattern 
...2121I121I12I.... 

• If the quasiholes are associated with the pat- 
tern . . . 2121il2l!l21 . . . , then after exchange the 

11' ° 

state will stay in same topological sector with 
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FIG. 11. Bratteli diagram of the fe = 3 Read-Rezayi state 
with two possible paths indicated. The dashed hne (red) 
corresponds to the sector label 12112112112112030212 and 
the dot-dashed line (blue) corresponds to the sector label 
302120302112112030. 



amplitude e^^ip~^ and transition with amplitude 
gi-7r(2-2a)(^-5 j^^^q ^ sector that has the ex- 
changed quasiholes associated with the pattern 
. . . 212l03030!21 .... 



These rules make it easy to visualize what is going on 
as a result of braiding in this non-Abelian state, as de- 
picted in Fig. |10| It remains to see which of our two 
solutions, if any, agrees with the representation of the 
braid group obtained from conformal block monodromies 
[7(Jj. To make contact between these various representa- 



tions, let us observe that our representation of topological 
sectors as patterns separated by domain walls is in nat- 
ural one-to-one correspondence with the representation 
given by paths meandering through a Bratteli diagram. 
Fig. 11 Here, the vertices of the diagram are associ- 
ated with the various ground-state patterns according to 
their "height" in the diagram, and the links represent the 
possible domain walls between them. A left-to-right path 
along the links of the diagram then represents an allowed 
sequence of patterns separated by domain walls, hence, 
a topological sector. The same diagrammatic labeling of 
sectors also naturally arises through fusion rules in the 
CFT analysis of the RR states. With this identification, 
it becomes easy to see that the above rules describing our 
solution are, for a = —1/5, in one-to-one correspondence 
with the "tensor representation" established in Ref. [701 
based on the analysis of conformal blocks. 

To make this point, we briefly review the latter. In 
the tensor representation given by Slingerland and Bais 
[70] , topological sectors, or paths through the Bratteli 
diagram, are represented by tensor products of vectors 

«A.,A, + i of the "domino" form l^Ai.Aa ® ^^Aa.Aa ® ^'A3,A4 ® 

■ ■ ■ <E) i'A„_i,A„- Here, represents the "height" of the 
i-th vertex in the path (Fig. 11), A^+i = A^ ± 1, and 



'I'Ai.Ai+i is a formal vector representing a link between 
two neighboring vertices at heights A,j, Ai_)_i, respectively. 
At general level k, A, takes on values 0, . . . , fc. 

In this tensor product basis, exchange of the quasiholes 
with indices i and i -I- 1 is represented by a matrix R^^i 
which acts only on the i-th and {i + l)-th factors flU\: 



Rk,i VK^^K^ + i ® UA, + l,A.+2 l^A,,Ai + l ® WA. + l,A,-|-2 

RkA WA,,Ai-l ® WA,-l.A.-2 = a^^A,,Ai-l ® VK,-lM--2 



(171) 



Rk,i WA,,Ai-|-l <8) WAi + l,A, = -aq ^WA.,A, + l WA,-|-l,Ai (A^ = 0) 
^fe,iWA,,Ai-l 8)'yA,-l,A, = -ag^^WA^.A,-! «) WA,-l,Ai [K = k) 



(172) 



Rk,iVKiM + l®Vf^. + 'i-M \ ^ 
Rk,i VA- (»VA^-l^Ai I 



LA, + 1J, 
1 VLA.+2J,LA,J, 
LAi + lJ, 



-aq 



,1 VLA.+2J,LA.J, 



aq 



LA. + IJ, 



«A.,A. + 1®«A. + 1,A. (0<A,<fc) 
VAiM-1 8) WAi-l,Ai / 

(173) 



where q = , a = qW^M+rj ^ m jg related to the filling sequences of patterns via paths in the Bratteli diagram, 
factor via v = 3/(3Af + 2), and "g-deformed integers" 
[mj q are defined as. 



To wit, our first rule is equivalent to Eq. (171). To see 



\m\q = 



12 — a 2 



(174) 



this, observe that the two domain walls defined by the 
ground-state sequence . . . 303021211121 . . . could b e rep- 
resented on the Bratteli diagram Fig. ITT] by Wq i (g) Wi.2 
or by V3^2 ^ V2,i, both of which follow the form of 
Eq. (171). A similar observation can be made about 
. . . 212l|l212|030 .... The phase picked up by the states 

([T73| reduce to the rules established in the beginning of in Eq. \l7l\ is a = , which is equivalent to the phase 



In our case k = 3 and M = 0, so q = e J' and a — eio . 



In this case, it is not difficult to check that Eqs. (171) 



this section, once tensor products are reinterpreted as in our first rule e™il+'^<^)^ where a 



-h here and in 

5 



33 



the following. Similarly, one can observe that the pat- 
tern in o ur second rule is represented by the vectors in 
Eq, 
whic 



172 1. The phase in that equation is — ^ = 



is equivalent to the phase in the second rule, 
. Finally, our t hird and fourth rules are together 
equivalent to Eq. ([l73|. The patterns . . . 21210303^21 . . . 



and . . . 212l!l2lll21 . 

I I 

wi,2 ® « 2,i o r as U2,3 ' 



. can be written as vi^ <S vo,i and 
) U3,2 and V2.1 ® Vi^2, which appear 
1 and Ai = 2, respectively. Up to 
a change in the order of the basis states, the matrix in 
Eq. ( 173 ) for either value of Aj gives the matrix elements 



in Eq. (173) for A^ 



stated in the third and fourth rules; this equivalence is 
shown here for A,; = 2: 



-aq 



"9 i3i; 

,1 ^/l4j;l2j, 

L3J, 



"9 L3J, 



aq 



L3J, 



3 If 



e 



:> If 
lif- 



gi7r(i+a)j^-l 



^i7r(|— 2a) 



T(l-2a)^-i 



(175) 



We hence see that one of our two solutions does indeed 
agree with the prediction based on conformal block mon- 
odromies, with the other one being closely related. 

Furthermore, it appears that the solutions we obtained 
form a true subset of the solutions that can be derived 
by imposing the relevant fusion rules, together with the 
axioms defining general anyon models (see, e.g., Refs. 
021 US [71]). If, in addition to the pentagon and hexagon 
equations, one imposes unitarity and modularity, these 
admit four solutions [HI [T^j. Two of these appear to 
be identical to ours, with the other two related to the 
former by complex conjugation. We observe that in our 
approach, there is no reason to expect that solutions au- 
tomatically come in complex conjugate pairs. This is 
so since the coherent state ansatz explicitly assumes a 
holo morphi c dependence on quasihole coordinates (see 
Sec. Ill A I, corresponding to a choice of sign for the 



magnetic field that renders trial wave functions for the 
RR state holomorphic (in both electron and quasihole 
coordinates). Our findings thus seem to imply that for 
the "missing" two solutions, one cannot construct holo- 
morphic trial wave functions that can be adiabatically 
deformed (through a continuous family of local Hamilto- 
nians) into the thin torus patterns we work with. 



VI. CONCLUSION 

In this work, we have further developed an approach 
to construct the braiding statistics associated to a given 
fractional quantum Hall state through adiabatic trans- 
port of quasiparticles. This approach is based on the 
notion of adiabatic continuity between FQH on the torus 
and simple product state — or "patterns" — in the thin 
torus limit. We have demonstrated that this notion, to- 
gether with a suitable coherent state ansatz for local- 
ized quasihole states, allows one to work out the result 
of adiabatic transport, using locality arguments and the 
information contained in the patterns. The latter in- 
clude properties under translation and the transforma- 
tion of sectors when particles are rearranged along cer- 
tain topologically nontrivial paths on the torus. The ap- 



proach also makes heavy use of modular invariance. We 
have presented a refined and unified treatment of simpler 
cases studied earlier in this formalism, and then moved 
on to demonstrate the applicability of the approach to 
the fc = 3 Read-Rezayi state. In all cases, we found 
results consistent with conformal block monodromies of 
trial wave functions. This is of particular interest in the 
RR case, where to our knowledge, there are no rigorous 
results yet that guarantee the agreement of these mon- 
odromies with the statistics defined through adiabatic 
transport. Our approach also has the benefit of giving 
rise to intuitive pictures representing the transformation 
of topological sectors during braiding, and allows for a 
logically independent, self-consistent derivation of non- 
Abelian statistics without heavy mathematical machin- 
ery. 

The findings presented in this work make us hopeful 
that the coherent state approach developed here, and ear- 
lier, quite generally yields a sufficiently constraining set 
of equations that allows one to construct the statistics 
of a given state of interest. This should be possible, at 
least for the large class of states that can be described by 
the recent paradigm of thin torus or dominance patterns. 
We have also tested this in detail on the gaffnian state 
[73| . This may be of particular interest since this state is 
associated with a nonunitary CFT, and has been argued 
to have gapless excitations in its bulk spectrum [TH [7S] . 
The general question of a well-defined notion of braiding 
statistics is therefore quite subtle. However, the question 
of whether or not our approach yields well defined statis- 
tics in this case, and whether they are consistent with 
conformal block monodromies, is well posed. We found 
the answer to be affirmative, and will present details of 
the calculation elsewhere. 
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There are two solutions to this equation: 



Appendix A: fc = 3 Read-Rezayi solution: Two 
quasiholes 



^12 = Ap2^n 
or = A^ii . 



(A5a) 
(A5b) 



We begin with Eqs. (141 1 and (142 1, the locality and 



unitarity conditions, respectively, and seek to constrain 



the ^2 ^, coefficients. As written, Eq. ( 142 1 does not pro- 
vide information about the overall phase of which is 
the overall phase relation between the two mutually dual 
bases. This phase is, a priori^ arbitrary. We have, how- 
ever, chosen a phase convention by defining the action 
of the antilinear operator r for both bases (in agreement 
with the phase relation chosen in Eq. (4|). The sym- 

which we can use 
obal pat h" alo ng x) 



We will now show that the first of the above equations 
never leads to consistent independent solutions, except in 
the special case ^13 = 0. To see this, we feed Eqs. (A5) 
back into Eq. (A2e), and find, respectively, that 



or ei3' = (P + P"')Cii'- 



(A6a) 
(A6b) 



metry under r gave rise to Eq. ( 136 ) 
together with Eq. ( 129[ ) (from the "g' 
to replace (C"*")^ in favor of rewriting Eq. (142) as 



We first utilize the above to study all cases with ^13 = 
0. This implies either — 0, or p e {i, —i}. The 
former leads to a contradiction in Eq. (A2c). It is then 



straightforward to show that for p e {i, — i}, the solutions 
of the system ( A2 1 produce the braid matrix 



^ 






A 











A 





e ' 2 = 13^3 . 




-iJ 
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(Al) 

We will ex pand this matrix product by plugging the form 
Eq. ( 132 ) for i^"'" derived from global path constraints. 




(A7) 



This gives four independent constraint equations for the 



with the upper (lower) sign corresponding to Eq. (A5a 
(Eq. (A5b|). Equation (A7) corresponds to a consistent 
solution to the constraint equations ( A2 1. However, when 





- AK12' - Ap'ei3' 


= 


(A2a) 


61 [Veil 


+ Ap^i2] + p£,13^33 


= 


(A2b) 




2Ap^n^i2+P^u 


= e^5 


(A2c) 




2^631^ ~ ^33^ 




(A2d) 



Eq. (A7) is generalized to an n-quasihole system using 

it is not difficult 



HIE 



Recall p = — exp \ in{s + ^)] and A'^ — 1. Similarly, we 
may use Eq. ( 136| i n the definition of the reduced braid 
matrix x, Eq. (140 1, writing x Expanding 
the latter again with Eq. ( 132 1 and comparing the result 



the locality arguments of Sec. 
to see that the resulting braid matrix violates the Yang- 
Baxter equation. While this might suffice to rule out 
this solution, we have emphasized in the beginning that 
our approach requires no a priori assumption that any 
aspect of quasiparticle exchange is topological. We will 
thus show more directly in App. |B]that Eq. (A2| leads 
to contradictions in the present framework when three 
quasiholes are considered. Since we can rule out the spe- 



to the locality constraint Eq. (141), we find two addi- 



tional independent constraint equations: 

(i+p')aia2-Apei3' = o 

£.31 Kii - ApCi2] + 6363 = 



cial solution leading to the upper sign in Eq. (A7|, this 
case has not been mentioned in the main text. 

We now proceed by exploring solutions with ^13 7^ 0. 



(A2e) 
(A2f) 



We first show that Eq. (A5a) does not lead to further 
independent solutions. To this end, we plug Eqs. (A5) 
first into Eq. (|A2f|), 



We will use the constraint equations ( A2 1 to solve for the 
unknown elements of x, the dots in Eq. ( 141 ), which can 
be found from the expansion of the product e"'"(e+)"^ 
be: 



to 



^13^33 

or 6363 



-{l~p')£3l£ll 
<1-P)£3l£ll, 



and similarly into Eq. (A2b) 



ei3 



P^£l3 



£11 + ^12 

P Cii + 4i2 - 

2^31^+^33 
A61 [ Vai + AC12] - Apei363 



Xll 
X22 
X33 
X23 



(A3a) 
(A3b) 
(A3c) 
X32 (A3d) 



C13C33 
or C13C33 



(P^ -p)C3iCii 

-(1-^)^31^11- 



(A8a) 
(A8b) 



(A9a) 
(A9b) 



While Eqs. (A8b) and (A9b) are identical, Eqs. (A8a 



We can break the solution of Eqs. ( A2 1 into two major 



sections, which are based on the two ways to satisfy the 





equation we obtain by combining Eqs. (A2e) and (A2a 
P ?ii 



and (A9a) turn out to be consistent with one another 
only in cases where both sides vanish on both equations. 
We have alr eady discus sed a ll cases with ^13 = 0. To 

we may thus focus on the 



and (A9a 



^12' - Aeiia2 - Ap2^ua2 



(A4) 



satisfy Eqs. (A8a 
case ^33 = 0. On the right-hand side, we can rule out 
^31 = because, with ^33 = 0, it contradicts Eq. ( A2d|. 
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We can similarly rule out 6i = because, with Eq. (A5) 
and Eq. (|A6]), it viol ates Eg. ( |A2c[) 
to solve both Eqs. (A8a) and ("Kd: 



The only other way 
I is to have p = ±1. 



In this case, however, both equations (A5) are identical. 
Thus, Eq. ( A5a[ ) does not produce independent valid so- 
lutions, except for p g {«, —i}, leading to the braid matrix 
Eq. ( A7 ) (upper sign) . As mentioned, the latter leads to 



inconsistencies in the case of three quasiholes. 



To find the solution to the constraints ( A2 ) that will be 



consistent with the three-quasihole case, we now discard 
Eq. ( A5a ) and proceed to work from Eq. ( A5b ) , and the 



equations (A6b), (A8b) derived from it. First, we plug 
Eqs. (A5b| and (A6b| into Eq. (A2c), which will give us 
an explicit form for ^n^: 



We need only plug into Eq. (Al3d) the square root of 
Eq. ( |A12 1 to write each element oi x in terms oi^^n^, for 
which we have the expression in Eq. (AlOl. We can also 



absorb the A factor in Eq. (A13d| into the ± induced 



by taking this square root. Thus we reach the following 
form of the braid matrix 



/p-i 

p{p + p-^-i) ±v/pTp^(i-p) 
\ ±VpTp^(i-p) p + p-^-i 

(A14) 



which was presented in the main text as Eq. ( 143 1 



(AlO) 



(l+p)2 ■ 

In particular ^ 0. From Eq. ( |A8b ) we thus obtain 



f 2^33 



(All) 



= (P + P-'K33' 



means of Eq. ( A2d ) then gives 



where we have used Eq. ( A6b ) . Elimination of ^33 by 



i. 



(P + P 



e 2 



31 - - /(i+p)2 
iP + p-')^ll' 



(A12) 



We can now revisit the unknown elements of x- We 
rewrite Eq. ( A3 ) using the equations we have developed 
above. 



p-\i+pr^^,' 



Xii 



(A13a) 
(Al3b) 

X33 (A13c) 
A(l -p)(l +p)2e3iai = X23 = X32 (Al3d) 



p{p + P ' -l)(l+p)'6i' -X22 

ip+p-' -m+pf^i^' 
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3 Read-Rezayi solution: Three 
quasiholes 



Here we will solve a system of equations for the ele- 
ments of the three-quasihole transition matrix elements 
S,a,a' and the resulting braid matrix. The procedure is 
the same as that employed for two quasiholes: Using 
various constraints on ^'^ already derived in the main 
text, we write out the matrix elements of the unitarity 
equa tion, ^++(C^^)^ =l4x4, and the locality constraint 
Eq. ( 144 1 . This gives a system for the remaining un- 



known elements of However, we must recall that the 
form in Eq. ( 144 ) was based, in part, on the two-quasihole 
braid matrix Eq. (143). In App. [Ajwe found one other 
"special" solution for this matrix, namely Eq. (A7) (up- 
per sign), that was not presented in the main text. Here 
we will consider this special solution also, giving rise to a 
modified version of Eq. ( 144 ) , and show that this solution 



leads to inconsistencies with three-quasihole braiding. 

Just as we did in App. |Aj we first use the (anti- 
linear) mirror symmetry to eliminate complex conjuga- 
tion from the definition of the reduced braid matrix, 
X = ^^^{S.^'^Y- This is achieved by using Eqs. ( |149[ ) 
and (157). The result is 



x = e^ 



f A^p 



-1 



(e 



++\T 



-1/ 



A- 



(Bl) 



1/ 



We will expand this matrix product using the constrained form of in Eq. ( 161 ), reproduced here 







63 


64 


A2ei3 


AVCu 


-A-V63 


-a-VCm 


63 


-A-VC13 


AV61 


-AKi4 


^41 


-A-VC41 


-AK4I 


C44 



(B2) 
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where we recall A* = 1, and p is defined in terms of 
the shift parameter s as before. We gain a system of 
constraint equ ations for the ^a,a'S by plugging Eq. (B2) 
into Eq. (Bl) and equating the product to one of the 



3 1^14^ + 1^44!^ = 1- By comparison with Eq. (B5e|, this 
implies 



(B6) 



following expressions for x that have been derived from 
locality and from consistency with the two-quasihole so- 
lution. For generic parameter p, we found that the latter 
must be of the form Eq. (144), which we reproduce here 



Once the ^a,a' are known, Eq. (Bl|) 
the following expressions for the un 



allows us to obtain 
fuown elements of x- 



Xii = X22 



(B7a) 



^loc 



v 



p{p + p ^ - 1) ±\/p~- 
±y/p+p^il-p) 



p-l(l-p) 

p + p^^ — 1 

(B3)' 

with blanks denoting zeros. The 2x2 block in the 
above was taken directly from the two-quasihole solu- 
tion, Eq. (A14), as explained in the main text. For 
p S however, we found an additional solution 

to the two-quasihole system of equations, leading to the 
form of the braid matrix Eq. ( A7 ) . Using this form and 



(2A-Veiiei3 



14 



X33 



-An^ + 2A^pCi3^-AVe 



e 



X44 
X34 = X43 



44 



(B7b) 
(B7c) 
(B7d) 



(^41 



the same reasoning that lead to Eq. (B3), for p S {z, —i} 
the reduced braid matrix must be of the form 



A- VCii + SA^Kis + Ap^i4U4 ] e-"^+'^ 

the 



If we subtract Eq. (B5b) (times A ) from Eq. (B5a 



^loc 



V 



TP 



(B4) 



resultant equation can be solved two ways. 

£.13 = -Ap^ii 
or = A^ii 



(B8a) 
(B8b) 



1/ 



where the lower sign is just a special case of Eq. (B3|, the equations 



We can quickly eliminate one of these possibilities by 
comparing to the equations from unitarity. When Eqs. 
(B8) are put into Eq. (B5g), they respectively produce 



but the upper sign corresponds to the "special" solution. 



We equate xioc Eq. (Bl|. We first focus on those 
matrix elements for which X]nr identically zero in all 



cases. By means of Eq. (B2|), this gives rise to the fol- 
lowing three equations: 



A(p - p')£is' + A- V^n' - ^P'£u = (B5a) 



or 



Ml 



141 



-3ieiir 

\£ii\'{p+p-'-i) 



(B9a) 
(B9b) 



the same equation for both cases: 



Whereas we can put Eqs. (B8I into Ec^. (B5d) and get 
in for both cases: 

ia4|' = -3|eii|' + l, (BIO) 



A-vei3' + (p'-p')ai63 



A-Vei4'=0 (B5b) 



41 



which clearly contradicts Eq. (B9a). Thus Eqs. (B8a 



P Cii+A (p-p )£i3 +A ^14^44 = 0. (B5c) and (|B9a[) are not true. Eliminating |^i4|^ from Eqs 



We will also use equations gained from enforcing the uni- 
tarit y of £~^~^. If we expand ^''"''(C'''^)^ = I4x4 using 
Eq. ( B2 ) we find the following independent equations: 



(B9b) and (BIO) gives us an expression for j^nj , which 
we can turn into an expression for ^^^^ with the inclusion 
of some phase 61. 



I6ir+2iei3r + i^i4r = 1 
2 



41 



31^41 

^11* - 2Ap^i3 



+ ie44r = 1 

Cl4*'?44 = 



A^asCii* +p'ai?i3* - A-vi63i'-A-ipiCi4i 



(B5d) 

(B5e) 

(B5f) 

^ = 
(B5g) 



Cii^ 



e'^^p 
JlTpf 



(Bll) 



Furthermore, putting Eq. ( B8b ) into either Eq. ( B5a 
Eq. (|B5b[) gives: 



or 



(B12) 



For convenience, we may also write the unitarity con- 
dition in the form ('C^''')^^'''^ = I4x4j yielding a simi- 
lar (and equivalent) set of equations, one of them being 



which fixes the phase between (^^^^ and ^n^. Together 
with Eq. ( B6 ) , this also implies 



e4i^=e^^^^A2(p + p-i-l)Cii^ 



(B13) 



37 



where we introduced another phase 62- Assu ming first 
that ^ 0, we define e'^^ = Ui/^u (cf. Eq. ([B6|), we 
may solve Eq. (|B5c|) for ^44: 



^44 = {2p-l)A^e'<'-^n- 



(B14) 



It is easy to see that the last equation also holds in cases 
where ^14 = C41 = 0.[71] With Eqs. (iBSbl), (IbT^I), pTsl), 



and (B14) we can rewrite the unknown elements of x in 
terms of £,11"^, and Eqs. (B7) become 



X44 
X34 



2ie2 



(B15a) 
(B15b) 
(B15c) 

(B15d) 



-i7rX-\-i7r 



Or, using Eq. (Bll), 



V 



■v) 



■V 



where we have defined e*^ = e*^^' 
quoted in the main text as Eq. 



- 1 e 

(B16) 

"^+". This is the result 
It is worth noting 
The 



(162) 



that once again, t he S p arameters have dropped out. 
derivation of Eq. (Bll I is valid irrespective of the value 
of p, since we did thus far not use the diagonal matrix 
elements of Xloc "^^ich may take on special values for 
p G We are now able to rule out p G {i,— i}, 

and thus the "sp ecial" solution obtained in App. |Aj For 
in this case, Eq. (B16) has non-zero off-diagonal matrix 
elements, whereas Eq. (B4) does not. This justifies Eqs. 
(1431 and (144) in the main text, which ignore the "spe- 



cial" solution. Requiring consistency between the non- 
zero matrix elements of Eqs. (144 1 and (162 1, which we 



have not done in this Appendix, finally provides informa- 
tion about the phase p, relating it to the golden mean. 
This short argument is presented in the main text. 
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